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ABSTRACT 

New CO interferometer data show the molecular gas in IR ultraluminous 
galaxies is in rotating nuclear disks or rings. The CO maps yield disk radii, 
kinematic major axes, rotation speeds, enclosed dynamical masses, and gas 
masses. The CO brightness temperatures, the double peaked CO line profiles, 
the limits on thermal continuum flux from dust, and the constraint that the 
gas mass must be less than the dynamical mass, all indicate the CO lines are 
subthermally excited and moderately opaque (r = 4 to 10). We fit kinematic 
models in which most of the CO flux comes from a moderate-density, warm, 
intercloud medium rather than self-gravitating clouds. Typical ring radii are 
300 to 800 parsecs. 

We derive gas masses not from a standard CO-to-mass ratio, but from a 
model of radiative transfer through subthermally excited CO in the molecular 
disks. This model yields gas masses of ~ 5 x 10^ Mq, ~ 5 times lower than 
the standard method, and a ratio M^i^/Lqq 0.8 M0(Kkms~-'^ pc^)~^. In the 
nuclear disks, we derive a ratio of gas to dynamical mass, Mgas/M^yn ~ 1/6 and 
a maximum ratio of gas to total mass surface density, /x//Xtot, of 1/3. 

For the galaxies VII Zw 31, Arp 193, and IRAS 10565+2448, the CO 
position- velocity diagrams provide good evidence for rotating molecular rings 
with a central gap. 

In addition to the rotating central rings or disks a new class of star 
formation region is identified which we call an extreme starburst. They have a 
characteristic size of only 100 pc. with about 10^ M0 of gas and an IR luminosity 
of ~ 3 X 10^^ Lq from recently formed OB stars. Four extreme starbursts are 
identified in the 3 closest galaxies in the sample including Arp 220, Arp 193 
and Mrk 273. They are the most prodigious star formation events in the local 
universe, each representing about 1000 times as many OB stars as 30 Doradus. 

In Mrk 231, the C0(2-l) velocity diagram along the line of nodes shows 
a l".2-diameter inner disk and a 3"-diameter outer disk. The narrow CO 
linewidth, the single-peak line profile, the equality of the major and minor axes, 
and the observed velocity gradients all imply the molecular disk is nearly face-on, 
yielding low optical and UV extinction to the AGN. Such a geometry means 
that the molecular disk cannot be heated by the AGN; the FIR luminosity of 
Mrk 231 is powered by a starburst, not the AGN. 

In Mrk 273, the CO (1-0) maps show long streamers of radius 5 kpc(7") 
with velocity gradients north-south, and a nuclear disk of radius 400 pc (0.6") 
with velocity gradients east- west. The nuclear disk contains a bright CO core of 
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radius 120 pc (0.2"). 

In Arp 220, the CO and 1.3 mm continuum maps show the two "nuclei" 
embedded in a central ring or disk at p. a. 50° and a fainter structure extending 
7" (3kpc) to the east, normal to the nuclear disk. Models of the CO and dust 
flux indicate the two K-hand sources contain high-density gas with n(H2) = 
2 X 10"^cm~^. There is no evidence that these sources really are the pre-merger 
nuclei. They are more likely to be compact, extreme star burst regions containing 
10^ Mq of dense molecular gas and new stars, but no old stars. Most of the HCN 
emission arises in the two "nuclei" . The luminosity to mass ratios for the CO 
sources in Arp 220 are compatible with the early phases of compact starbursts. 
There is a large mass of molecular gas currently forming stars with plenty of 
ionizing photons, and no obvious AGN. The entire bolometric luminosity of 
Arp 220 comes from starbursts, not an AGN. 

The CO maps show the gas in ultraluminous IR galaxies is in extended 
disks that cannot intercept all the power of central AGNs, even if they exist. 
We conclude that in ultraluminous IR galaxies — even in Mrk 231 that hosts 
a quasar — the far IR luminosity is powered by extreme starbursts in the 
molecular rings or disks, not by dust-enshrouded quasars. 

Subject headings: galaxies: nuclei — galaxies: interstellar matter — galaxies: 
starburst — galaxies: ISM: dust, extinction — galaxies: individual (Arp 193, 
Arp 220, Mrk 231, Mrk 273, VII Zw 31) — radio lines: galaxies 
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1. INTRODUCTION 

The infrared ultraluminous galaxies are recent galaxy mergers in which much of the 
gas in the former spiral disks has fallen into the center. The gas re-establishes rotational 
support at a radius of a few hundred parsecs. This is the same size scale as the Narrow Line 
Region in active galactic nuclei, and the rotating nuclear disks in some elliptical galaxies. 
During its infall from 5kpc to 0.5 kpc, the high-density gas forms many massive stars, 
as in the early starbursts that formed the bulge stars in protogalaxies at high redshifts. 
The ultraluminous galaxies we observe in the local universe are re-enacting this primeval 
process, with a large part of the central mass in gas ready to turn into stars — and ready 
to fall into a black hole. 

The large nuclear concentration of molecular gas in ultraluminous galaxies has been 
detected in the millimeter lines of CO by many groups during the past decade. Scaling to 
the signal strengths from Milky Way molecular clouds, however, soon led to a paradox for 
many of the sources — the estimated gas mass was equal to or larger than the dynamical 
mass indicated by the linewidths. For Arp 220, for example, Scoville et al. (1991) found 
that nearly all of the mass in the central few hundred parsecs was in the form of molecular 
gas. To resolve this dilemma, we showed that in the extreme environment in the central 
600 pc of ultraluminous galaxies, much of the CO luminosity must come from an intercloud 
medium that fills the whole volume, rather than from clouds bound by self gravity, and 
therefore the CO luminosity traces the geometric mean of the gas mass and the dynamical 
mass, rather than just the gas mass (Downes, Solomon, & Radford 1993). This allows the 
nuclear disk gas to be overluminous in CO, relative to self-gravitating clouds in galactic 
spiral arms. The important linewidth for this calculation is the dynamical linewidth of the 
galaxy, not that of individual self-gravitating GMCs. The physical basis for estimating gas 
mass from the CO emission is therefore different from that in Milky Way GMCs. The Milky 
Way conversion factor is relevant for an ensemble of GMCs in an ordinary spiral galaxy, 
but not for the center of an ultraluminous galaxy. 

In a subsequent survey paper on 37 ultraluminous galaxies observed with the IRAM 
30 m telescope (Solomon et al. 1997), we derived the molecular gas mass by several methods. 
The minimum estimates of the gas mass, for optically thin CO, did indeed reduce the gas 
mass to much lower fractions of the dynamical mass. The true gas mass must be between 
the optically thin CO estimate and the dynamical mass, both of which are less than the gas 
mass derived with the Milky Way conversion factor for self-gravitating molecular clouds. 

The significance of the lower molecular gas mass is twofold. First, it makes the gas 
mass smaller than the dynamical mass, as it must be. Second, it means the total molecular 
gas mass in ultraluminous galaxies is similar to, and not greater than, the molecular mass 
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in the disk of a gas-rich spiral. It is thus not necessary to force into the nuclear region all 
the gas in the galaxy — H2 and H I. The nuclear starbursts can be fed by pre-existing 
molecular gas that falls into the central few hundred parsecs from an original radius of a 
few kiloparsecs, and not from the outer H I disk. 

These ideas show that the kinematics of the nuclear gas is critical to understanding 
both the evolution of ultraluminous galaxies and the molecular line formation. The CO 
emissivity per unit gas mass surface density is not constant, but varies with position in the 
nuclear disk. A more correct treatment requires radiative transfer modeling, with kinematic 
data from millimeter interferometers. 

To carry out this next step in our study of molecular gas disks at the centers of IR 
ultraluminous galaxies, we made high-resolution CO maps with the IRAM interferometer. 
We chose sources from our CO survey, concentrating on the nearest ultraluminous galaxies 
in the northern sky, for which we already had accurate measurements of the total flux, 
linewidth, and radial velocity. The new observations give accurate positions, line and 
continuum fluxes at the interferometer resolution, images, directions of the kinematic 
major axes, and radii of the nuclear disks. Q These radii and the velocity extrema give 
the enclosed dynamical mass. The velocity data allow modeling of the kinematics and 
the CO brightness temperature versus radius. For most sources, the CO data indicate 
the gas kinetic temperature is higher than in typical Galactic molecular clouds, but the 
CO is subthermally excited, turbulent, and warm, quite unlike the molecular clouds in 
galaxy-scale disks. At the positions of maximum CO intensity in the disks, the CO lines are 
moderately opaque, with optical depth r(l-O) ^ 4 to 10, but in other parts of the disk, the 
CO lines can be optically thin. The result is that the molecular gas mass is about a factor 
of five lower than it would be if the CO emission came from self-gravitating clouds. The 
rest of the dynamical mass is accounted for by young stars that provide the luminosity and 
older stars that were already in the central bulge before the merger. 

We argue that the bright spots in the ultraluminous galaxies are compact extreme 
starburst regions that occur in the dense gas traced by HCN and CS emission. They 
produce most of the HCN luminosity, and 1mm dust continuum while the more diffuse gas 
in the nuclear disks dominates the CO luminosity. These compact extreme starburst regions 
which we identify in Arp 193, Arp 220 and Mrk 273 are the most prodigious star formation 
events in the local universe. There are undoubtedly many more too compact to be resolved. 

The plan of this paper is as follows. Section 2 describes the observations and 
summarizes the parameters obtained by direct measurement. Section 3 introduces our 



^In this paper, distances are for Hq = 75kms ^ Mpc ^ and qq = 0.5. 
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model of a turbulent rotating disk and discusses the derivations of the rotation curve, disk 
radius, disk height, turbulent velocity, gas mass, and dynamical mass obtained by fitting 
the model disk to the data. Sections 4 through 8 present the results on the galaxies with 
CO (2-1) data, obtained with sub-arcsecond resolution. These galaxies are VII Zw 31, 
Mrk 231, Arp 193, Mrk 273, and Arp 220. Section 9 gives the results on five additional 
ultraluminous galaxies mapped in CO(l-O) only. Section 10 discusses the source sizes, gas 
masses, the ratio of gas mass to dynamical mass, and the mass of young stars needed to 
power the starburst. Section 11 contains our conclusions. 

2. OBSERVATIONS AND RESULTS 

CO(l-O) and (2-1) were observed with the IRAM interferometer on Plateau de Bure, 
France (Guilloteau et al. 1992). Three to four configurations of the four 15 m antennas 
gave baselines from 24 to 410 m for VII ZwSl, Mrk 231, Mrk 273, and Arp 220, and 24 to 
288 m for the other sources. The SIS mixers had receiver temperatures of 60 - 80 K, and 
operated single sideband at 3 mm and double sideband at 1.3 mm with system temperatures 
of 200 and 400 K, respectively. The spectral correlator covered a 500 MHz band at 2.5 MHz 
resolution, giving velocity ranges of 1300 km s~^ at 7kms~''" resolution at 3 mm, and 
700 kms^^ at 3.4kms~^ resolution at 1.3mm. For analysis, we smoothed the data to 20 
and 40kms~^. Amplitude and phase were calibrated relative to quasars and the fiux scale 
was adopted from same-epoch measurements of quasars and planets with the interferometer 
and the IRAM 30 m telescope. The beamwidths (FWHM) were 1" to 3" at CO(l-O) and 
0".5 to 1" at C0(2-l). Table 1 summarizes the observing parameters. 

Table 2 lists CO source positions and integrated CO line fiuxes from the interferometer. 
We also list the CO line fiuxes from the 30 m telescope (Radford, Solomon, & Downes 
1991a; Solomon et al. 1997). Table 3 lists sizes from Gaussian fits to the CO visibilities. 
Sources observed in CO(l-O) and (2-1) had the same sizes in both lines. The visibility fits 
are independent of the synthesized beam or the CLEAN algorithm, and there is no need 
to deconvolve an apparent size on a map. The phase calibration can broaden a source due 
to the baseline error, varying as 27rX/D times the angle between the source and its phase 
calibrator (A = wavelength, D = baseline). For sources and calibrators in this paper, the 
broadening is < 0.1 beamwidth. After correcting for atmospheric decorrelation, we obtained 
sizes < 0".2 for the quasar calibrator sources. 



3. DISK MODELS 
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3.1. Geometry of Nuclecir Disks from Kinematic Data 



CO surveys of our Galaxy show most of the gas in the inner R < 5kpc is in ~ 6000 
giant molecular clouds (GMCs) of diameter ~ 50 pc and H2 density 150 cm~^. In an 
ultraluminous IR galaxy merger, these clouds fall into the central R ~ 500 pc, forming a 
disk of height ~ 50 pc. Because the GMCs' density is too low to stabihze them against 
tidal shear at the merger's center, and because the new volume is ten times smaller than 
the original volume of all the GMCs in the galactic disk, the GMCs loose their identity and 
blend into a continuous medium with a mean H2 density ~ 10^ cm~^. We therefore modeled 
the nuclear disks as a continuous medium rather than an ensemble of individual clouds. 

We used the CO spectral data to estimate how much of the gas was in a high-density, 
inner ring or disk and how much in a lower-density, outer disk. To simulate the observed 

spectra and position-velocity diagrams, we modified a model for rotating disks developed 
by A. Dutrey (see Dutrey, Guilloteau, & Simon 1994). This model has an inner disk of 
high-density gas between radii -Rmin and Ri = Rmm + and an outer disk of lower-density 
gas between radii Rmin and i?max- The H2 density between Rmin and i?max is 



where the inner disk is the Gaussian of width AR (FWHM), and the outer disk is the power 
law. We fit the outer disks with constant density {a — 0). Our beams cannot distinguish 

between a ring and a filled disk, so for most sources, we took /?min = 0. Because of the 
radiative transfer through the rotating, inclined disk, the predicted CO map looks like a 
ring, even if there is gas in to i? = 0. For the galaxies VII Zw 31, Arp 193, and IRAS 
10565H-2448, the position-velocity diagrams are better fit with ring models, with Rmin — 
i?o, rather than filled disk models. Table 4 gives the source geometry, derived as follows: 

Line of nodes: The visibility fits in individual spectral channels gave position offsets 
vs. velocity, which showed the direction of the kinematic major axis. We checked the angle 
of the line of nodes from channel maps and isovclocity contour maps. Errors are ±10°. 

Rotation curve turnover radius Rq: In the data, the position- velocity diagrams along 
the kinematic major axis often have two peaks near the source center. In our models, twin 
peaks occur near the points where the rotation curve turns over and becomes flat. The 
observed offsets of the two intensity peaks on the major axis were tried as first guesses for 
the rotation curve turnover radii Rq. Model fits converged close to these values. In angular 
units, the uncertainties are ~ ±0".l. In our models, Rq is the projected radius with the 
largest gas column density at the same hne-of-sight velocity, and the highest flux. For most 
of the sources, the sohd angle inside this zone is too small to affect the maps, even if the 
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gas is really in a filled disk rather than a ring. 

Inner disk half-intensity radius, Ri = Rmin + Ai?; In our models, most of the CO 
flux comes from the high-density inner disk. The extent of the central double peak in the 
position-velocity diagrams constrained the radial width AR, which we varied until the 
model position- velocity diagrams matched the data. The uncertainty in AR is ~ ±0".2 

The outer radius, Rmax, of the low-density disk was taken to be the observed CO 
maximum extent on the line of nodes. 

Disk thickness H: Along the z-axis perpendicular to the equatorial plane, the gas 
volume density n{z) near z = was approximated by a gaussian with full width to 
half-maximum H. For a disk with a fiat rotation curve, the disk's thickness can be 
estimated with the Mestel formula (Mestel 1963; see also Binney & Tremaine 1987): 



H{R) = 1.4 a{R)(^ 



'^■^W^e^U.ILT' , (2) 



V, 
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Pi V 2RJ 



where a is the 1-D velocity dispersion (assumed independent of z) and pi{R) is the total 
mass density (gas plus stars) in an equivalent sphere with a flat rotation curve. This 
equivalent density, in MqPc"^, is related to the rotation velocity, in kms~^, through the 
usual integral for dynamical mass, pi{R) — 18.5(Kot/-R)^- We dropped the second-order 
term and estimated the disk thickness as 



H ^ 1.4 a{R) 
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(3) 



where H and R are in pc, Kot is in kms~^, Pgas is in MqPC~^. Because the velocity 
dispersion decreases with distance from the center, our model disk thickness tends to a 
constant value. In reality, the molecular disks are warped and twisted, and their shape 
depends on the merger's infall history. 



3.2. Rotation Curves 

We combined our density models with rotation curves that reproduced the observed 
channel maps, position-velocity diagrams, spectral profiles, interferometer visibilities, and 
fine intensities. We took the rotation velocity to be 

V.ot{R) = Vo (-j , (4) 

with P — 1 ioT R < Rq and P — ior Rq < R < -Rmax, that is, the curve rises from the 
center and flattens after Rq. Most of the derived values of Rq are smaller than our beam, 
so more elaborate rotation curves are not justified for now. 
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For the inclination, i, we started from the major/minor axis ratio of the integrated 
CO source, with i = cos~^ (minor/major), and then varied the inchnation to match the 
position- velocity data and to ensure the derived gas mass was less than the dynamical 
mass. The nuclear rings or disks may not be circular; they may be elliptical, analogous to 
X2 orbits at the center of a barred galaxy. The inclination i listed in Table 5 should be 
interpreted in the sense sin i = sin / cos 0, where / is the true inclination of the disk relative 
to face-on, and is the azimuth of the major axis of the eUipse, if there is one. With our 
current beams, we cannot tell the difference between circular and elliptical orbits. 

For the rotation velocity, Vq, we first guessed the apparent speed, Vosini, to be half the 
velocity range of the twin CO peaks on the major axis. The best fits were close to these 
values, giving 200 to 300kms~^ (corrected for inclination) on the flat part of the curve 
(Table 5). This is an equivalent circular velocity if the orbits are elliptical. 

The turbulent velocity, AV , of the gas was taken from the model fit, not from the 
observed line profiles, which include both turbulent broadening and the rotational velocity 
gradient in the beam. Rotation alone cannot explain the observed profiles. We need a local, 
turbulent line broadening as well. The models actually constrain this parameter rather well. 
We assumed the local fine broadening had the form 



where AV is the local fine halfwidth to the 1/e level {AV = 0.6 x FWHM = lAa, where 
a is the r.m.s. velocity dispersion along the line of sight). The radiative transfer model 
included the local linewidth at each point and the rotational velocity versus radius. We 
iterated the model fits until the local turbulence and the rotational velocity gradient — 
convolved with our beam — matched the observed position-velocity diagrams and the 
observed line profiles. Turbulent velocities AV that fit the data were 30 to 140kms~^, with 
some of the position- velocity diagrams indicating lower turbulence at greater radii from the 
nucleus. At the rotation curve turnover point Rq, the values of AV correspond to local 
FWHM hnewidths of 70 to 230kms~^, so the nuclear disks are highly turbulent. 



After deriving sizes and velocities from the data, we adjusted our models to match 
the observed CO fine intensities. As a first guess, we set the gas kinetic temperature 
equal to the dust temperature derived from blackbody fits to the IRAS fluxes. At the 
gas densities in our models, collisions cannot raise the CO excitation temperature to 
the gas kinetic temperature, nor populate the CO levels thermally. Instead, there are 
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3.3. CO line peirameters 
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different excitation temperatures for each transition. The hne brightness temperatures, 

excitation temperatures, and gas kinetic temperatures are in the sense T{, < Tex < T\^m, with 
Tcx(2 — 1) < Tex(l — 0). We therefore assumed the CO level populations were determined 
by the local gas kinetic temperature, density, and velocity gradient, and we calculated, with 
an escape probability program, CO column densities per unit velocity width for a local 
velocity gradient of lkms~^pc~^ and a [CO/H2] abundance of 8 x 10~^, both typical of 
Milky Way molecular clouds. For these values, the escape probability method yields the 
following non-LTE excitation temperatures T^^: 

rex(l - 0) = 0.186 n°-622 ^0^240 ^g) 

rex(2 - 1) = 0.117 r°n'' , (7) 

where n is the H2 number density in cm~^, and Tkin is the gas kinetic temperature in K. 
These formulae are valid for 50 < Tkin < 150 K and 100 < n< 3000 cm'^. 

We then calculated the radiative transfer on the hnes of sight through the disk (eqs. 7 
& 8 of Dutrey et al. 1994), convolved the computed map with our beam, and adjusted the 
parameters to match the observed brightness temperatures. Highest CO opacity occurs at 
the rotation curve turnover radius, Rq, which yields the largest column density at the same 
line-of-sight velocity. Table 6 gives the best-fit model CO(l-O) excitation temperatures 
and gas densities, and Tables 7 and 8 give the model CO(l-O) and (2-1) line opacities and 
Rayleigh- Jeans brightness temperatures. We also ran models for higher density gas, with 
excitation temperatures closer to the kinetic temperatures, but these models gave too much 
CO and dust flux, and did not reproduce the position velocity diagrams. 

At the rotation curve turnover radius Rq, our models yield CO(l-O) opacities of 2 
to 8. These mean opacities apply to smoothly distributed gas. The model of subthermal 
excitation yields the observed CO(2-1)/(1-0) ratios of 0.6 to 1.0 (Radford et al. 1991a). In 
the lower-density, outer disks detected in these sources, the volume-averaged densities in 
the model are too low to yield detectable CO hnes. For these outer disks, our model must 
be corrected with an area filling factor for molecular clouds, as in CO maps of disks of 
normal galaxies. For the mean values for the outer disks in Table 6, we adopted a volume 
filling factor of 0.1, an area filling factor of 0.3, and local gas densities 10 times higher than 
the volume-averaged values listed in the Table. 

3.4. Gas Mass and Dynamical Mass 

We used the CO parameters to estimate the gas mass and the dynamical mass vs. 
radius (Table 9). After we found a density law that reproduced the observed brightness 
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temperatures, we integrated over the source to get the gas mass. These masses are for a 
[CO/H2] abundance of 8 x 10~^, as in molecular clouds in our Galaxy. For the sources 
observed here, it is difficult to lower the model CO abundance below the Milky Way value. 
Keeping the same CO line intensity at a lower CO abundance would force us to make the 
gas mass bigger than the dynamical mass, invalidating the model. 

The dynamical masses were taken to be RV^^^/G, where V^ot is the rotation speed 
listed in Table 5. Most of the sources have fiat rotation curves, so the dynamical mass 
increases linearly with radius. Table 9 lists the ratio of gas mass to dynamical mass in the 
inner, high-density disk alone (to the radius Ri in Table 4) and in the inner and outer disks 
together (to the radius -Rmax in Table 4). Relative to our assumed Hubble constant, the gas 
mass scales as Hq^^ the dynamical mass as Hq^ and their ratio as Hq^. 

We also list in Table 9 the maximum ratio of gas to total (gas + stars) surface density 
in our models. We calculated the ratio of gas to total mass surface density, n/ntot, in 15 
bins in radius, and found maximum values of ~l/3 to 1/4, midway between the rotation 
curve turnover radius Rq and the half-intensity radius Ri. The models in this paper are for 
the distributed gas, but there is also denser gas in star-forming cores that give rise to HCN 
and CS lines. Most of the CO flux comes from the distributed medium, but our gas masses 
may have to be corrected upward to allow for the cores, depending on future interferometer 
results in the dense-gas tracer lines. We return to this point in the Discussion section. 



3.5. Dust Continuum Flux 

We used the model gas mass to predict the dust flux, from 

= (1 + z)K{ur)MDl'B{ur, To) , (8) 

where S is the dust continuum flux density, i^obs and Vr are the observed and rest-frame 
frequencies, M is the gas mass, Dl is the luminosity distance, B is the Planck function, 
and Td is the dust temperature. We took the dust temperatures from our fits to the IRAS 

fluxes (Solomon ct al. 1997) — the same as the gas temperatures in Table 6 — and we used 
a dust mass absorption coefficient k{v,.) = 0.1 x z/", where z/^ is in THz, k, is in cm^gm^^ of 
interstellar matter, and we took the index n = 1.5. At 230 GHz, this absorption coefficient 
is K = 0.011 cm^ gm~^ of interstellar matter, or, for a gas-to-dust mass ratio of 100, Kd = 
1.1 cm^ gm~^ of dust, as in previous estimates for dense molecular clouds (see Kriigel & 
Siebenmorgen 1994, their Fig. 12). The predicted thermal fluxes from dust are then 

5(i/obs) = 6.4 X 10-^(1 + z)M Df Ta 1/,'+" , (9) 
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where flux density S is in Jy, gas mass M is in Mq, luminosity distance Dl is in Mpc, dust 
temperature is in K, and rest frequency i^r is in THz. The predicted thermal dust fluxes 
at 2.6 and 1.3mm are compared with the observed fluxes in Table 10. 

Four sources, VII Zw 31, Mrk 231, Mrk 272, and Arp 220, were observed in both the 
CO(l-O) and (2-1) hues, and on the longest basehnes of 410m. Because we have better 
angular resolution for these sources, we discuss them in more detail than our other sources. 

4. VII Zwicky 31 

This source, catalogued by Zwicky (1971), was identifled by Fairclough (1986) as an 
ultraluminous IR galaxy. Its optical surface brightness proflles resemble those of elliptical 
galaxies and its optical hne ratios flt a starburst rather than an AGN (Djorgovski, de 
Carvalho, & Thompson 1990). Sage & Solomon (1987) found this source had one of the 
highest known CO luminosities and deduced that the molecular gas was a large part of the 
dynamical mass. Scoville et al. (1989) showed the CO source had a radius < 3.8 kpc. 

The CO maps of the galaxy VII Zw 31 are the best evidence in our sample for a 
rotating ring. Models that best fit the data are those with a minimum ring radius i^min 
in eq.(l) equal to Rq = 290 pc. Filled-disk models, with i?mm = 0, give poorer fits, with 
much less contrast. In 40kms~^-wide channel maps, the CO peak migrates from south to 
north with increasing velocity (Fig. 1). The source is resolved east-west, perpendicular 
to the kinematic major axis. This indicates the source is inclined, but more face-on than 
edge-on. Our modeling gives a good fit to the data for a ring inclined at 20° to face-on. 
Figure 2 shows the maps of integrated intensity and isovelocity contours across the source 
in both CO lines. The total CO(l-O) flux measured with the interferometer is equal to 
the flux obtained at the 30m telescope (Radford et al. 1991a). At C0(2~l), the source is 
partially resolved out, and the flux from the interferometer is only half the single-dish flux. 
In channels off the CO line, the continuum flux is < 2mJy at 109 GHz and < lOmJy at 
221 GHz. This limit on any thermal dust flux is consistent with the gas mass of 1.1 x 10^° 
deduced from the CO luminosity, for a dust temperature < 50 K. 

The CO velocity contours in Fig. 2 indicate the north-south velocity gradient. In the 
position- velocity diagrams (Fig. 3) there is a 200 km s^^ velocity shift over 4" north-south. 
The model in the Tables reproduces well the observed diagrams, the CO profiles across 
the source, the channel maps, and the observed line intensities. The CO spectra (Fig. 4) 
have twin peaks separated by ±70kms~^. The blue and redshifted CO peaks are 
respectively 0".65 south and north of the source centroid. The hne profiles are remarkably 
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symmetric along the north-south kinematic major axis. The CO hne of nodes differs from 
the northwest-southeast dust lane found by Djorgovski et al. (1990) on a 10" scale by 
subtracting a model from their optical images. 

5. Markarian 231 

The Seyfert I galaxy Markarian 231 has an IR luminosity of 3.3 x 10^^ (e.g., Sanders 
et al. 1987). Most of the power is emitted at the center of a major galaxy merger which has 
tidal tails extending over 75kpc (Cutri, Rieke, & Lebofsky 1984; Hutchings & Neff 1987; 
Sanders et al. 1987). Optical spectroscopy indicates stars formed in the past 1 Gyr over 
a widespread region (10-15 kpc). In this large extranuclear region the Ha emission comes 
from shocks rather than stellar photoionization (e.g., Hamilton & Keel 1987; Lipari, Colina, 
& Macchetto 1994). The best previous CO study is that by Bryant & Scoville (1996), who 
found an east-west velocity gradient in the CO source. Our higher-sensitivity results have a 
factor of two better resolution than their study, and not only confirm the velocity gradient, 
but also distinguish the inner and outer disks in the position- velocity diagrams. 

5.1. CO in Mrk 231 's molecular disk 

In contrast to the large optical extent, the CO source is very compact, on a scale twenty 
times smaller in radius. In 40kms~^ channel maps in CO (1-0) and (2-1), the line peak 
shifts from west to east with increasing velocity (Fig. 5). The gradient is clearly shown on 
the CO isovelocity maps (Fig. 6). The C0(2-l) intensity contours (Fig. 6) are symmetric 
and only slightly broader than the beam. The CO (2-1) position- velocity diagram along 
the line of nodes (Fig. 7) shows an inner nuclear disk of diameter 1".2 (radius 460 pc) and 
a 3" outer disk with a lower velocity gradient. The asymmetric C0(2-l) line profiles 0".4 
from the center of the galaxy (Fig. 8) are a characteristic signature of a rotating nuclear 
disk. The maximum linewidth in both CO lines is 190kms~^ FWHM, which is narrow for 
ultraluminous galaxies. 

The CO data imply the molecular disk is face-on, for three reasons. 1 ) The CO major 
and minor axes are nearly equal, indicating an inclination i < 20°. To reproduce the 
single-peaked CO profile, narrow CO linewidth, and observed velocity gradients, our model 
disks must be within 20° of face-on. For i > 20° to face-on, the gas mass would exceed the 
dynamical mass. 3) The CO does not absorb the nuclear continuum source at any velocity. 
This agrees with the CO disk being face-on, leaving the Seyfert I nucleus unobscured. For 
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1 — 10° to face-on, the true rotation velocity would be 345kms~^, the same as we deduce 
for Arp 220. If the disk were even more face-on, the observed velocity gradient would imply 
a rotation speed > 400kms~^, higher than in all the other galaxies in our sample, which 
seems unlikely. We adopt i = 10° and Vrot = 345kms~^, although the kinematic data alone 
would also allow a rotation velocity of 250kms~^. 

In CO, we do not detect the gas emitting in the 1667 MHz OH megamaser with 
760 km s""*^ linewidth (Staveley-Smith et al. 1987). nor do wc detect anything at the IR 
source 3". 5 south of the nucleus, to a limit of 5mJy in 40kms^-'^ channels (sec also Bryant & 
Scoville 1996). Armus et al. (1994) interpreted the southern IR source as the nucleus of the 
merger partner, but recent HST images show this feature to be a dense arc of star-forming 
knots (Surace et al. 1998). In any case, it is very weak in both CO and the mm continuum. 

In the center of Mrk 231, the AGN is detected by centimeter wavelength VLBI as a 
variable, nonthermal radio continuum source with a size < 1 pc (Preuss & Fosbury 1983; 
Neff & Ulvestad 1988; Lonsdale, Smith, & Lonsdale 1993; Taylor et al. 1994). New VLBA 
images at 1.4 GHz resolve this nonthermal source into a nuclear core with pc-scale lobes. 
Carilli, Wrobel & Ulvestad (1998) subtracted their VLBA map from a VLA map and found 
an 0".4 nonthermal nuclear disk source that emits half the nuclear flux, the rest coming 
from the core and the pc-scale lobes. 

At millimeter wavelengths, outside the CO line channels, we detect an unresolved 
(size < 0".l) nonthermal, mm continuum source at the nucleus, which is probably the 
nonthermal VLBA core. The mm continuum source and the CO centroid both coincide 
with the 8.44 GHz continuum source (Condon et al. 1991). 

5.2. Mrk 231's Luminosity: 2/3 Starburst, 1/3 AGN 

The CO kinematic data for the central 1 kpc of Mrk 231 and our model of a face-on 
disk show the far-IR luminosity comes from a starburst, not a black hole accretion disk. 
We assume the central UV-to-FIR continuum has three components: the visible and UV 
continuum is blackbody flux from the AGN accretion disk, the 2fim peak is blackbody flux 
from a dusty AGN torus, and the far-IR peak is blackbody flux from the molecular disk. 
The extinction to the accretion disk is low. Ay = 2 mag, because both the dusty torus and 
the larger-scale molecular disk are face-on. 

The accretion disk's luminosity may be derived from the UV and visible fluxes and 
the extinction to the Seyfert I nucleus. The UV fluxes measured by lUE and HST are 

2 X 10"^^ erg s"^ cm'^ A"^ (Schmidt & Miller 1985; Hutchings & Neff 1987; Smith et al. 
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1995). A blackbody with this flux, corrected for a UV extinction of 5 mag (Smith et al. 
1995) and a distance of 170 Mpc, has a UV luminosity of 4 x 10^^ Lq. For an accretion disk 
temperature of 20000 to 30000 K, as in the models of Malkan & Sargent (1982) and Sanders 
et al. (1989), the blackbody peaks at 2000 to 1000 A, and has a radius of 1 x 10~^pc. This 
result is consistent with the optical continuum, which has Ay = 2 mag and an intrinsic 
magnitude My = —25.1 (luminosity 5 x 10"*^^ Lq; Boksenberg et al. 1977). 

The dusty torus' luminosity may be derived from the near-IR continuum. As in other 
Seyfert I galaxies, the 2fim bump is well fit by a 1470 K blackbody, close to the dust 
sublimation point (Kobayashi et al. 1993). The flux at the 2/xm bump, 9.3 x 10~^^ erg s~^ 
cm^^/xm^^, corrected for an IR extinction of 0.2 to 0.6 mag (Krabbe et al. 1997), implies 
a blackbody luminosity of (3 to 4) xlO^^L©, and a radius of 0.18 pc. The compact 10 
to 25 urn source in Mrk 231 (Matthews et al. 1987; Keto et al. 1992; Miles et al. 1996) 
arises farther out in the dusty torus, where the dust temperature is ~200K, at a radius of 
~10pc from the AGN. This dusty torus is heated by the AGN. If its thickness is equal to 
its radius, then the torus absorbs half the AGN's power, so the AGN's total luminosity is 6 
to 9 X 10^^ Lq, the same as estimated from the optical and UV flux. 

The molecular disk's FIR luminosity may be derived from the IRAS fluxes, and is 2 x 10^^ Lq 
(e.g., Solomon et al. 1997). The molecular disk has a dust optical depth close to unity 
at 100 /xm and causes the deep lOfim silicate absorption (Allen 1976; Rieke 1976; Roche, 
Aitken, & Whitmore 1983). The blackbody radius calculated from the FIR luminosity is 
200 pc, half the measured size of the molecular disk. The molecular disk receives some 
heat from the dusty torus. With the thickness-to-radius ratio in our CO model, it would 
intercept half of the power, or 2 x 10^^ Lq of re-processed radiation from the AGN. Since 
this is only 10% of the FIR luminosity, most of the FIR luminosity of Mrk 231 must come 
from the starburst in the molecular ring or disk. A similar model was used by Rigopoulou, 
Lawrence, & Rowan-Robinson (1996) to flt the FIR continuum spectrum. 

Our GO model flt gives a gas mass of 3.1 x 10^ Mq — equal to the mass of all the 
molecular gas in the Milky Way — in a high-density disk between R — 76 and 850 pc. 
There is another 0.9 x 10^ Mq of gas in a lower density outer disk extending to a radius of 
1.7kpc. The high-density molecular disk is the region observed in H2 at 2/xm by Krabbe et 
al. (1997). The brightest H2 line is 200 pc from the nucleus, in the molecular disk, so the 
low 2 /im extinction on the line of sight to the nucleus itself does not apply to the H2 flux. 

Our conclusion that most of Mrk 231's FIR luminosity comes from a starburst can be 
checked for consistency (Table 11). We assume the new stars in the burst, Af„cw*- yield 
Lfir/Miicw* ~ 500 Lq Mq^^, a luminosity ratio that can be attained in a fast starburst 
(Leitherer & Heckman 1995). Dividing this ratio into the FIR luminosity yields Mnew*, 



-16- 



from which we then calculate M^id-k — M^y^ — Mg^g — M^^w*, where Moid* is the mass of 
old stars in the miclcar bulge before the merger. For the region R < 460 pc in Mrk 231, 
the rotation velocity inferred from our CO data yields a dynamical mass Mdyn(< 460 pc) 
= 12.7 X 10^ M© , and our CO model yields Mgas = 1-8 x 10^ M©. The consistency check 
then yields Moid* = 7.1 x 10'' M0, which is about the same as the mass of old stars in 
a similar radius in the Milky Way (e.g., Oort 1977). Because of our assumed L/M, the 
derived stellar masses in the starburst and the bulge are uncertain by at least a factor of 
two. Our point is simply that with a starburst powering the FIR luminosity, all the gas and 
stellar masses are quite plausible. In summary, we think Mrk 231 's black hole accretion disk 
emits ~ 1 X 10^^ L©, or ~ 30% of the total bolometric luminosity. Most of the luminosity, 
~ 2 x 10^^ Lq, comes from the starburst in the molecular ring or disk. 



6. Arp 193 

The ultraluminous merger Arp 193 shows two tidal tails in the visible and the near 
IR (Smith et al. 1995, 1996). The IR images show a single, elongated nucleus, but 
cm-radio continuum maps show two sources separated by 1". Strong lines of Br7 and 
H2 are seen in K band (Goldader et al. 1995). The ncar-IR continuum is dominated by 
young supergiants formed in the merger-induced starburst. Previous CO line ratio studies 
indicated sub-thermal CO excitation and moderate densities (Radford et al. 1991). Our 
detection of HCN in Arp 193 indicated, however, that large amounts of dense gas were also 
present (Solomon, Downes, & Radford 1992). 

The C0(2-l) emission is nearly unresolved in individual 20kms^^ channels, with 
temperatures up to 8 K (Fig. 9). The CO intensity and isovelocity contour maps (Fig. 10) 
reveal a rotating disk with a line of nodes at the same orientation as the optical isophotes in 
the center of the merger. The position- velocity diagrams and CO spectra (Figs. 11 and 12) 
along the line of nodes have the characteristic signature of a rotating ring, with two peaks 
near the rotation curve turnover radius, as in VII Zw 31 and IRAS 10565-1-24. Models 
that best fit the data are those with a minimum ring radius i?min in eq.(l) equal to Rq — 
220 pc. As in VII Zw 31, fiUed-disk models, with -Rmin = 0, give poorer fits, with much less 
contrast. The ccntroid of the CO source coincides with the main radio continuum peak on 
the 8.4 GHz map by Condon et al. (1991). The strongest CO peak is 1" to the southeast of 
the ccntroid, and peaks at — 120kms~^ (relative to 225.282 GHz). It is nearly unresolved on 
our maps, coincides with a secondary radio continuum peak. This hot, compact, southeast 
CO core in Arp 193 is responsible for the large difference in the CO(2-1)/(1-0) ratio in the 
blue and redshifted sides of the fine profiles we observed at the 30 m telescope (Radford et 
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al. 1991) . This region is hotter and denser than the disk and has properties similar to a 
huge molecular cloud core. 

This compact southeast core is one of several sources we identify in this study as 
Extreme Starhurst Regions (see Section 8 on Arp 220 for a detailed discussion). Using the 
observed CO luminosity, linewidth, radius, and temperature, we can estimate the gas mass, 
dynamical mass and far IR luminosity of the core. The results, summarized in Tabic 12, 
indicate an object with a gas mass ~ 6 x lO^M©, a luminosity of ~ 2 x 10^^ Lq, and about 
~ 1 X 10^ Mq of newly formed stars. 

7. Markarian 273 

The galaxy Markarian 273 has at least one Seyfert 2 nucleus with optical linewidths 
of 700kms~^ (e.g., Sargent 1972; Koski 1978). There is evidence for two Seyfert 2 nuclei 
(Asatrian, Petrosian, & Borngen 1990), and two nuclei are seen on IR images (Armus et 
al. 1992; Majewski et al. 1993; Knapen et al. 1997). Centimeter radio maps show a 0".3 
northwest peak coinciding with the CO, and a weaker 0".2 southeast peak with no CO or IR 
counterpart (Ulvestad & Wilson 1984; Schmelz, Baan, & Haschick 1988; Sopp & Alexander 
1991; Condon et al. 1991). There is also an 18cm VLBI source of 16mJy (Lonsdale et 
al. 1993; Smith, Lonsdale, & Lonsdale 1998a). The best previous CO study is by Yun & 
Scoville (1995), who observed the source with a 2" beam and found a component extended 
north-south as well as an unresolved nuclear component. Our higher-resolution CO maps 
also show the extended structure, but now resolve the nuclear source into a nuclear disk 
oriented east- west, and a very compact core embedded in the nuclear disk. 

7.1. Extended molecular gas and nuclear disk in Mrk 273 

The extended gas is best traced in the CO(l-O) channel maps, blueshifted to the south 
and redshifted to the north (Fig. 13a). The CO at -1-250 to -1-350 km s~^ runs in an arc 
from 3" to 7" (2 to 5kpc) north of the nucleus (Fig. 14). This extended gas also appears in 
the maps of CO(l-O) integrated intensity and isovelocity contours (Fig. 15a). The CO(l-O) 
lines are weaker and narrower in the northern arc than in the nuclear disk (Fig. 16). The 
arc has 4% of the CO(l-O) flux of Mrk 273, so if gas mass is proportional to CO flux, then 
the northern arc has a gas mass of 1 x 10^ Mq. Figure 17a is the CO position- velocity cut 
in declination, through the nucleus. The narrow-line emission extends 7" (5 kpc) south at 
Okms"^ and 7" north at -1-300 km s"^. It is tempting to speculate that the streamers are 
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bringing molecular gas into the center. The north-south extent of the streamers on our CO 
map agrees with that on the map by Yun & Scoville (1995). While those authors deduced 
a line of nodes at p. a. 30°, our isovclocity contour maps show there are two perpendicular 
kinematic systems. In the extended gas the velocity gradient is north-south, and in the 
nuclear disk it is east- west. 

The 2" nuclear disk with its compact core are best seen in the CO (2-1) channel maps, 
made with a 0".6 beam (Fig. 13b). The nuclear disk kinematics are shown in the CO (2-1) 
map (Fig. 15b), where the isovelocity contours trace the east-west velocity gradient of the 
nuclear disk over a 1 kpc diameter. Figure 17b is the CO position-velocity diagram of 
the nuclear disk in right ascension. The broad nuclear disk line has a full width to half 
maximum of 380 km s"^, and its velocity centroid changes from —300 to -1-200 km s~^ over 
0".5 from west to east. This is a projected velocity gradient of 1.5 kms~^ pc~^, the same 
as deduced by Schmelz, Baan, & Haschick (1988) from H I absorption. Table 11 lists the 
stellar and gas masses that would be needed in a star burst in the Mrk 273 molecular disk if 
the ivprn/Mnew* ratio is the same as we adopted in our starburst model of Mrk 231. 

7.2. An Extreme Starburst Region in Mrk 273 

The high-resolution C0(2-l) maps show a remarkable molecular-hne source in the 
Mrk 273 nuclear disk — a bright, 0.35" x < 0.2" CO core (Fig. 15b), that resembles the 
west nucleus of Arp 220. This is the most luminous extreme starburst region in our sample 
of 10 galaxies. It has an IR luminosity of about 6 x 10^^ Lq , generated in a region with a 
radius of only 120 pc and a current molecular mass of 1 x 10^ Mq (Table 12). To put this 
in perspective, the entire molecular core has a radius about 5 times that of an IR luminous 
Milky Way CMC (for example, W51) but with about 3,000 times the molecular mass and 
^ 10^ times the IR luminosity from OB stars. 

This core has a broad CO line with a zero-intensity width of 1060 km s^^, the same as 
the OH megamascr (Stavclcy-Smith et al. 1987). It coincides with the northwest extended 
continuum peak at 8.44 GHz, which has a size 0".32 x 0".18 (Condon et al. 1991), the 
same as the compact CO. The radio spectrum is nonthermal, with a spectral index of —0.6. 
The thermal dust emission dominates the spectrum above 350 GHz (Chini et al. 1989; 
Rigopoulou et al. 1996). At the CO core, we detect continuum fluxes of ll±2mJy at 
111 GHz and 8 ± 2mJy at 225 GHz. Extrapolation of the synchrotron and thermal dust 
spectra indicates the dust contributes 50% of the flux at 1.3 mm and < 10% at 3 mm. The 
extended nonthermal continuum emission coincident with a high mass of dust and gas 
leaves little doubt that this region is powered by star formation. 
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8. Arp 220 

The center of the ultraluminous galaxy Arp 220 has two radio continuum and two 
IR sources 1".0 apart that are interpreted as merger nuclei. These are not exactly the 
same objects in the radio and the IR, because both of the east and west IR nuclei are four 
times larger (0".8) than the radio sources (0".2) (Condon et al. 1991; Graham et al. 1990; 
Majewski et al. 1993; Miles et al. 1996; Scoville et al. 1998). The X-band continuum is 
starlight associated with the nuclei, and is best seen on the HST NICMOS images (Scoville 
et al. 1998). These images resolve the eastern IR nucleus into a strong NE and weaker SE 
component, the latter of which coincides with the eastern radio continuum component and 
the eastern OH masers (Diamond et al. 1989). The two radio sources are extended and 
nonthermal, and are produced by supernovae in the most active star-forming regions. The 
1.3 and 1.6 /im [Fe II] lines in the west source indicate iron evaporated from dust in shocks 
(Armus et al. 1995a; van der Werf & Israel 1998). The best previous CO study is that 
by Scoville, Yun, & Bryant (1997), who derived the kinematics of the nuclear disk with a 
resolution of 1". 



8.1. CO and dust emission from the two nuclei of Arp 220 

Our C0(2-l) and 1.3mm continuum maps (Fig. 18a) show two compact sources 
embedded in more extended emission. The continuum fluxes at 1.3 mm are due to dust, 
because they are well above the extrapolated radio synchrotron spectra (Fig. 20). The 
different appearance of the dust continuum and the CO maps is a temperature effect. At 
1.3 mm the dust is optically thin and its flux varies as the column density and the dust 
temperature, while the CO flux is partly opaque, with flux varying with temperatiire 
at the surface, and partly optically thin, with flux varying inversely with CO excitation 
temperature. The Arp 220 molecular disk contributes strongly to the CO maps, while the 
warmer compact peaks dominate the dust continuum maps. From the 1.3 mm dust flux 
alone, and for dust at 100 K and solar metallicity, we derive gas masses for Arp 220-west 
and -east of 0.6 and 1.1 x lO^^M^. 

In Arp 220-west, the 0".3 (100 pc) mm continuum dust source coincides within 0".l 
of the cm-continuum west peak, and the i^-band west peak (Scoville et al. 1998). In the 
CO (2-1) channel maps (Fig. 19b), Arp 220-west appears as a strong source in the range 
—310 to -1-210 km s~^, peaking at — llOkms"^ (czisr = 5340 km s~^), with an observed 
C0(2-l) brightness temperature of 27 K, of which 10 K is from the more extended Arp 220 
disk. The positive- velocity emission near Arp 220-west is centered at -|-80kms'"^, and 
appears in both the C0(2-l) channel maps and in a position-velocity cut at p. a. 300° 
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(Fig. 22b). This positive- velocity emission has about 25% of the CO integrated intensity of 
the negative-velocity emission of Arp 220-west. 

The dust continuum is centered between the negative and positive velocity emission, 
suggesting that Arp 220-west is a composite structure, with a kinematic axis nearly 
perpendicular to that of Arp 220's main disk. At the negative velocity west core, the 
emission extends asymmetrically from —50 to — 400 kms"-*^ (Fig- 22b). At the positive 
velocity core, it runs from +50 to +350 kms^^. These extreme, asymmetric line profiles, 
confined to a small region, suggest possible radial fiows along a bar-like structure, or 
immense high-velocity molecular outflows similar to those observed in high mass star 
formation regions. Molecular outflows are also the likely sites of the strong H2 lines (Sturm 
et al. 1996) . 

In Arp 220-east, the CO covers a wide velocity range. The negative-velocity gas at 
—230 to — 30kms~^ is separated from Arp 220-west by 0."85 at p. a. 110° — the same angle 
as the cm-radio continuum sources, the 18 cm OH masers, and the X-band SE peak (Scoville 
et al. 1998). The center velocity of this gas is czisr — 5330 km s~^, close to that of the 
OH megamasers and H2CO emission at the cm-radio east source (Baan & Haschick 1995; 
Lonsdale et al. 1998). There is also more extended, positive- velocity gas (0".9 FWHM) at 
130 to 290 kms-^ It is 1".3 from Arp 220-west at p. a. 85° — the same orientation as the 
iiT-band NE source (Scoville et al. 1998). This positive- velocity gas in Arp 220-cast, at 
~ +200kms~"^ {cz\sr = 5650 kms^^), has the same velocity offset from Arp 220-west as the 
ionized gas seen in Br7 and Pa/? (Larkin et al. 1995), and the weak OH masers at 1612 MHz 
(Baan & Haschick 1987). Its observed maximum C0(2-l) brightness temperature is 19 K, 
of which 4K is from the extended Arp 220 disk. The 1.3 mm continuum from dust in 
Arp 220-east has a diameter of 0".6 (205 pc) and is displaced by 0".3 from the cm-radio east 
source. It falls in between the X-band NE and SE peaks on the images by Scoville et al. 
(1998), and in between the positive and negative- velocity CO. The kinematic data suggest 
that Arp 220-cast has a steady progression of velocity between the CO peaks at 5330 km s^^ 
and 5650 km s^^, along the same position angle (p. a. 50°), and in the same sense of rotation 
as the main, larger-scale molecular disk of Arp 220 (Fig. 22c). 



8.2. CO in the Arp 220 molecular disk 

The molecular disk is best seen as an extended source in the CO(l-O) maps (Fig. 18b). 

The CO (2-1) maps also show this same extended source, when the compact peaks are 
subtracted. The disk has CO halfwidths of 2".0 x 1".6 at p. a. 50°, with a line of nodes 
at this same position angle, blueshifted in the southwest and redshifted in the northeast. 
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as seen in the CO spectra along the major axis (Fig. 21). Our model fits yield a rotation 
curve turnover radius of 200 pc, a disk outer radius of 480 pc. and a disk rotation velocity of 
330 km s^^ on the fiat part of the rotation curve. The CO emission centroid of the molecular 
disk is ~ 0".3 east of the western nucleus. These parameters reproduce the position- velocity 
diagrams along the kinematic major axis (Fig. 22), and are similar to the values derived by 
Scoville et al. (1997) by a different algorithm. On a larger scale, CO can be traced in the 
8" diameter outer disk at the same position angle. The inner and outer disks both contain 
shocked gas, detected in the near-IR vibrational lines of H2 (van der Werf 1996). 

The CO(l-O) maps also show an eastern streamer extending 7" (3kpc) perpendicular 
to the disk, in the range +45 to +245kms~^, curving in to the east nucleus with increasing 
velocity offset (Fig. 19a). This east streamer is most intense at 85 to 165kms~^ (Fig. 18b). 
It may be material that is still falling into the center. Surprisingly, this CO streamer 
appears on the HST F-band image by Shaya et al. (1994). The CO inner and outer disks 
coincide with the prominent optical dust lane at 50° (Fig. 23a), while the CO east streamer 
coincides with the perpendicular dust lane in the optical image (Fig. 23b). 

Our model fits to the CO position- velocity diagrams of the Arp 220 disk (Fig. 22) yield 
a much higher dynamical mass than estimates from the IR (e.g., Doyon et al. 1994; Shier, 

Ricke, & Rieke 1994, 1996; Larkin et al. 1995). From the masses in Table 9, we estimate 
M(< 350 pc) = 8.8 X 10^ Mq, and M(< 600 pc) = 1.5 x 10^° Mq, three to four times higher 
than the IR estimates. The millimeter CO data and the 2.3 fim CO bandhead absorption 
data yield different masses because the two IR nuclei do not lie on the line of nodes traced 
by the CO maps and because the near IR starlight is obscured and does not trace the full 
extent of the molecular disk. 

The gas mass derived from both the dust and the CO implies a column density of 
1 X lO^^cm"^ and Ay ~ 1000 mag through the disk, consistent with the dust having an 
opacity of unity at 180 /xm (Emerson et al. 1984; Scoville et al. 1991). The high opacity in 
the far IR may partly explain why the C"^ 158 //m line is weak in Arp 220 (Fischer et al. 
1998.). The ratios of fine structure lines observed by the ISO sateUite yield an equivalent 
screen Ay ~ 45 mag (Genzel et al. 1998). As is well known, a screen model gives only 
an extreme lower limit. The corresponding extinction when the emitting gas and dust 
are completely mixed is Ay ~ 1000 mag, the same as we deduce from the column density 
derived from the CO maps. 

Figure 24 shows our model of the molecular disk and the two "nuclei" . The CO disk 

is inclined 40° from face-on, Arp 220- west has a radius of 68 pc, and Arp 220-east has a 
radius of 110 pc. The disk thickness is 90 pc, which means the path to the near surface of 
the east and west K-band sources is to 20 pc, which is why the two nuclei are visible at 
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all at K band. That is, although the visible extinction through the entire molecular disk 
is 1000 mag. the visible extinction to the near surfaces of the east and west sources is only 
about 50 mag, because of the shorter path. Our diagram is similar to that of Scoville et 
al. (1997; their Fig. 7), except our presentation is a view from the pole of the disk. In our 
model, the near side is north, the far side is south. This agrees with the optical colors, 
which are bluer on the north side of the dust lane (Shaya et al. 1994). Since the CO east 
nucleus is south of the major axis, it is on the far side. The CO disk or ring includes both 
nuclei but the nuclei are oriented east- west, and hence the difference in their line of sight 
velocities does not give the full rotation speed of the disk. 

Our data and model have some important consequences for the interpretation of the 
Arp 220 nuclear sources. 

1) The rotation curve of the CO disk indicates a dynamical mass of 12 X 10^ Mq interior 
to 480 pc, which corresponds to the central bulge mass of a large spiral like the Milky Way. 

2) The gas mass in each of the two extreme starburst "nuclei" is only 6 x 10® Mq. Their 
individual luminosities are ~ 3 x 10^^ Lq. (About half of the Arp 220 FIR luminosity comes 
from the molecular disk, not the two nuclei). With an L/M ratio of 1000, corresponding 
to a super-starburst in its initial phase, the mass of new stars in each nucleus would be 
only 3 X 10® Mq, sufficient to explain all the K-band continuum luminosity. The velocity 
dispersion of the CO in each of the two nuclei imphes a dynamical mass of ~ 1 x 10^ Mq, 
about the same as our estimate of the sum of the gas and new stars. 

3) Hence, there is no room left over for old stars in the two "nuclei" — they cannot 
be the relicts of the old nuclei of the pre-merger galaxies. Furthermore, there is no 
observational evidence — radio, infrared, or optical — that they contain old stars. 

4) In any case, the masses of the two "nuclei" are negligible in comparison with the 
mass that controls the motion of the molecular disk. Furthermore, the two "nuclei" of 
Arp 220 have radial velocities indicating that they take part in the general disk rotation, 
i.e., they share the general rotation in the potential of the old bulge, and are dominated by 
its gravity, not their own. 

5) In our interpretation, the two "nuclei" of Arp 220 are not the pre-merger nuclei at 
all. We think they are just large ensembles of pre-existing GMCs that have been strongly 
compressed on their infall into the old bulge potential as a result of the merger. The masses 
of these large, compressed gas structures can reach 10^ Mq, and the strong compression, 
in about one orbital period, leads to the extreme starbursts, with luminosities of a few 
times 10^^ Lq, from each compact 100-pc structure. These extreme starbursts power the 
ultr aluminous galaxies. We are not seeing the birth of quasars. 
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8.3. Arguments for a steirburst 

What produces Arp 220's luminosity? We think it's stars, for three reasons: there is no 
obvious AGN, there are enough ionizing photons for a starburst, and there is enough dense 
molecular gas to make stars. 

1. Arp 220 has no AGN lines. The mid-IR lines of [O IV] and [Ne V] are not detected 
(Sturm et al. 1996; Lutz et al. 1996). The Bra width of 1300 km s^^ reported by Depoy, 
Becklin, & Geballe (1987) as evidence for an AGN was not confirmed in later measurements 
(Goldader et al. 1995; Larkin et al. 1995). 

2. Arp 220 has no AGN in the radio continuum. New VLBA maps (Smith et al. 1998b) 
indicate all the VLBI point sources are young radio supernovae in a dense medium, with 
radio luminosities comparable to SN1986J in NGC891. 

3. The radio nuclei are extended starburst regions. At centimeter wavelengths, their 
diameters are: west, 0".21 x 0".14 (72 x 48 pc) and east, 0".32 x 0".19 (110 x 64 pc) (Condon 
et al. 1991). The source extent, the low brightness at 5 GHz, and the high FIR-to-radio flux 
ratio all show the nonthermal radio continuum is starburst-dominated, not AGN-dominated 
(Sopp & Alexander 1991; Condon et al. 1991; Baan & Haschick 1995). 

4. Arp 220 is not a "warm" ultraluminous galaxy. IR ultraluminous galaxies with an AGN, 
like Mrk 231, Mrk 1014, and 08572+3915, have warm mid-IR colors {f 2b I ho > 0.2) and a 
bright, symmetric nucleus resembling a reddened QSO (Surace et al. 1998). Unhke these 
sources, Arp 220 has strong far-IR and weaker mid-IR flux (/25//60 — 0.08). It has been 
argued that a QSO is hidden in Arp 220, but the silicate absorption depth at 10 //m in 
Arp 220 (e.g.. Smith, Aitken, & Roche 1989; Dudley & Wynn-WiUiams 1997) has been 
greatly overestimated, due to the strong PAH lines directly adjacent to the silicate feature 
(Genzel et al. 1998), so a deeply embedded AGN is not required. 

5. The diffuse OH megamasers do not indicate black holes. A new VLBI study (Lonsdale 
et al. 1998) shows Arp 220 has two OH megamaser components in each nucleus, one diffuse 
and the other compact. The diffuse megamaser component arises in the extended starburst 
medium, and can be explained by IR pumping via photons absorbed in the 35 and 53 /im 
lines of OH. The observed 9- Jy depth of the 35 fim. OH line can only be explained by 
absorption against an extended source (Skinner et al. 1997). The com,pact OH megamasers, 
probably coUisionally pumped shocks, are currently the only evidence for possible AGNs in 
Arp 220 (Lonsdale et al. 1998), but it is hard to understand why such AGNs would have 
no associated radio continuum core-jet sources. 
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8.4. The ratio of bolometric flux to ionizing flux resembles that in Sgr B2. 

From the observed Br7 flux, Armus et al. (1995b) and Shier, Rieke, & Rieke (1996) 
estimated that < 10% of Arp 220's luminosity came from a starburst. Armus et al. noted, 
however, if they had underestimated the K band extinction, then a starburst could power 
the source. The higher near-IR extinction is indeed required by the high column densities 
found from the CO maps. Prom the radio free-free continuum, Scoville et al. (1991; 1997) 
claimed there were not enough ionizing photons for a starburst to power Arp 220. We 
think this expectation was wrong, because galactic star- forming regions like W49, W51, 
and Sgr B2 have FIR luminosities greatly exceeding the Lyman luminosities derived from 
their radio fluxes. In compact HII regions like W3(0H), or in our galactic center, the FIR 
excess is a factor of ~ 20 (e.g., Zylka et al. 1995). 

In H II regions with densities > 10^ cm^^, most of the Lyman photons heat dust rather 
than ionizing the gas (Jennings 1975; Panagia 1977; Fazio 1978; Mezger 1985), so the 
Lyman photon rate derived from the radio flux is only a lower limit. We may estimate this 
lower limit to Arp 220's Lyman continuum photon flux from the continuum at 113 GHz, 
which includes 15mJy from the synchrotron spectra and ISmJy from dust (Fig. 20). This 
leaves 13mJy as optically thin free-free flux, and implies a Lyman continuum photon rate 
> 1 X 10^^ s"^. Similar lower limits come from the H92q; radio recombination line (Zhao et 
al. 1996), and the low-frequency turnover of the radio spectra of the two nuclei (Fig. 20), 
due to free-free absorption in ionized gas with an emission measure of 10^ cm~^ pc and a 
mass of 5 X WMq (Sopp & Alexander 1991). 

For comparison, in the galactic giant H II region and molecular cloud Sgr B2, the 
FIR power is 10 times the ionizing luminosity derived from the radio (Gatley et al. 1978). 

Sgr B2 resembles Arp 220's molecular gas in column density, strong molecular lines, and 
far-IR opacity, but no one has ever claimed Sgr B2's FIR excess proves it is powered by 
a black hole. Arp 220 and Sgr B2 both have the same ratio of FIR flux (in W m~^) to 
radio free-free flux (in Jy), and hence the same FIR excess. So our lower limit on Arp 220's 
ionizing flux, scaled by the Sgr B2 template, yields exactly Arp 220's FIR luminosity of 
1.2 X 10^^ Lq, and is indeed compatible with a starburst. Similar conclusions are reached 
by Genzel et al. (1998) from ISO hue data on the ionizing flux. 



8.5. Arp 220 has plenty of gEis for a steirburst. 

Our models yield gas masses of 0.6 x 10^ M0 for Arp 220-west, 1 x 10^ M0 for 
Arp 220-east, and 3 x 10^ for the disk to R2 = 1.2 kpc (Table 9). The total gas mass 
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is 5 X 10^ Mq, about 6 times lower than previous estimates with a standard ratio of gas 
mass to CO luminosity. Our value agrees with the gas mass of 3.5 x 10^ Mq obtained by 
Sturm et al. (1996) from the H2 lines at 6.9 and 17/xm, although this latter value depends 
sensitively on the assumed temperature. Our method differs from that of Scoville, Yun, & 
Bryant (1997), who assumed all the dynamical mass is in gas, rather than stars, thereby 
forcing their model CO disk to be only 16 pc thick at R =250 pc, and to have a mean 
density of 1.5 x lO^cm"^. Because we do not assume all the mass is gas, we derive a disk 
thickness of 80 pc and a lower mean density than Scoville et al. Our higher- resolution data 
indicate the gas density is highest at the two nuclei, not at the center of the molecular disk 
as in the model of Scoville et al., and our value for the total gas mass is therefore a factor 
of two lower than their estimate. 

Even the lower molecular gas mass, however, is about the same as in the entire disk of 
a gas-rich spiral! It is enough gas for a huge nuclear starburst. In the Arp 220 inner disk, 
the gas surface density is 100 times the peak value in the Milky Way 4-kpc molecular ring 
and more than 20 times the surface density in the inner few hundred pc of our Galaxy. To 
estimate the stellar mass in the nuclear bulge, we adopted a global luminosity to starburst 
mass ratio in the Arp 220 disk of SOOLqMq^^, based on starburst L/M ratios from the 
models of Leitherer & Heckman (1995). This gives a estimate of the mass in newly-formed 
stars from the current merger-induced starbursts. Subtracting this mass and the gas mass 
from the dynamical mass then gives an rough estimate of the mass in old, pre-merger stars. 
The resulting bulge mass is 8 x 10^ M0 within R — 480 pc (Table 11), about the same as 
in a large galaxy like our own (e.g., Oort 1977). The starburst interpretation thus implies 
that old stars are a significant part of the dynamical mass. 



8.6. Arp 220-west has 25% of the total luminosity. 

The CO lines, the 1.3 mm dust flux, and the cm-radio continuum all suggest the 
currently most active starburst is Arp 220-west. Its radius of 68 pc and dust temperature 
of 75 to 80 K imply a blackbody luminosity of 3 x 10^^ L©, or about 25% of Arp 220's 
total IR luminosity of 1.3 x 10^^ Lq, as expected from the R^T^ ratio of the west and disk 
sources. This estimate also agrees with the fraction of the mid-lR flux absorbed by the 
35 ;um OH line at 5320 kms~^ (Skinner et al. 1997), the velocity of Arp 220 west. If the 
Lyman continuum scales as FIR luminosity, then Arp 220-west has 1000 times the ionizing 
flux of 30 Doradus, which is 10^^ photons s"^ (Kennicut & Chu 1994), from 2400 OB 
stars in a region of comparable size (Parker 1993). Because Arp 220-west is small and the 
dynamical time scale is short, its starburst must be younger (age 5 x 10^ yr) than that in 



-26- 



the larger Arp 220 disk. The initial phase of a compact starburst can be highly luminous, 

reaching L/Af ratios of 1000 - 3000 LqMq"^ (see models in Leitherer & Heckman 1995). 
Our estimate of the total mass (from the CO lincwidth) and the gas mass (from the CO 
luminosity), suggests that in Arp 220- west, 50 to 60% of the original gas mass has already 
turned into the new stars in the current, extreme starburst (Table 12). 

9. SOURCES OBSERVED IN CO(l-O) ONLY 

We observed these sources with lower resolution and thus our kinematic models are 
less precise. All these sources show the same phenomena as the previous group of sources, 
however, with velocity gradients and line profiles that suggest rotating rings or disks. 

9.1. 00057+4021 

The galaxy IRAS 00057+4021 has not been well studied, despite its high luminosity of 

LiR = 4 X 10^^ Lq, and its OH megamaser (Kazes, Mirabel, & Combes 1988). The i?-band 
image by Armus, Heckman, & Miley (1987) shows a 60" (48kpc) southeast-northwest disk 
with a tidal tail stretching a further 45" (36kpc) southeast. In contrast to the optical light, 
the CO source is quite compact. Most of the molecular gas is in a l".lx < 0".6 core source. 
The CO has wider lines (140 km s^^ FWHM) and is more intense on the southeast side of 
the fine of nodes than on the northwest (Fig. 25). The continuum flux from the CO source 
is <10mJy at 110 GHz. This hmit on the thermal dust flux is consistent with the gas mass 
of 9 X 10^ Mq deduced from our model flts to the CO. 

In the position-velocity diagram along the line of nodes at p. a. 135° the CO line has a 
velocity range of 300kms~^ over a diameter of 3" (Fig. 25). The CO fines are blue-shifted 
in the northwest and redshifted in the southeast. In strips parallel to the kinematic major 
axis the CO spectra show the double-peak behavior of a rotating disk or ring. To the 
southeast, the CO profile is steep toward the red, and gently sloping toward the blue. To 
the northwest, the profile is reversed. The nuclear velocity gradient in CO has the same 
major axis as the larger-scale optical isophotes. The position-velocity data and the spectra 
both indicate another, lower- intensity (0.5 K) disk with a broad linewidth (200kms~^) in 
the same region as the brighter, high-density disk. 
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9.2. 02483+4302 

The galaxy IRAS 02483+4302 is a merger with a tidal tail extending 90" to the west. 
The merger has two nuclei separated by 3". 8 east-west. In optical lines, nucleus A (west) 
has a Seyfert 2 spectrum, nucleus B (east) has a LINER spectrum. The optically more 
intense nucleus A appears to belong to an elliptical galaxy plowing through the disk of a 
former spiral containing nucleus B (KoUatschny et al. 1991). The optical continuum of 
nucleus B comes partly from hot stars (Womble et al. 1990). 

The CO (Table 2) coincides with nucleus B, supporting the idea that this was originally 
a gas-rich spiral. Our CO position also agrees within 0.5" of the ISmJy source seen at 
5 GHz (Crawford et al. 1996). At 109 GHz, the continuum flux from the CO source is 
< 10 mJy. This limit on the thermal dust flux is consistent with the gas mass deduced from 
the CO. In velocity channels 40kms^^ wide, the CO source is very compact. Figure 26 
shows the maps of CO integrated intensity and isovelocity contours across the source. The 
CO has a north-south velocity gradient, with a kinematic major axis perpendicular to the 
optical tidal tail extending west from the merger nuclei. 

In strips parallel to the kinematic major axis the CO spectra show the symmetric 
behavior of a rotating disk or ring. In the south, the CO profile is steep toward the blue and 
gently sloping toward the red. In the north, the profile is reversed. The position-velocity 
diagram along the line of nodes shows a gradient of 200kms~^ over 0".9 in declination 
(Fig. 26). The narrow CO linewidth and the low optical extinction to the nucleus both 
suggest the molecular disk around nucleus B is face-on. 

The CO source, at ^ = 0.05144, {czisr = 15420 km s"^) and the quasar Q0248+430 at 
z = 1.311 are separated on the sky by 16". 9 (the galaxy-quasar separation is incorrectly 
hsted as 3".5 by Burbidge (1996) and by Hoyle & Burbidge (1996)). The tidal tail that 
extends > 80 kpc from the two merger nuclei crosses the quasar at a projected distance of 
15.4 kpc from the CO source. At this position, low-density atomic gas in the tidal tail is seen 
in absorption against the quasar in the Na I D and Ca II H and K lines a,t z — 0.0515 and 
0.0523, with linewidths < 150kms~^ (Womble et al. 1990). There are two other absorption 
systems, in Mg II, at z = 0.394 and 0.451, from the halos of other, more distant galaxies 
on the line of sight to the quasar (Womble et al. 1990; Sargent & Steidel 1990; Borgeest et 
al. 1991). After correcting for the primary beam, we derived the quasar's continuum flux 
to be 190 mJy at 110 GHz in June-October 1994. Toward the quasar, our data at 20kms~^ 
resolution do not show any CO(l-O) absorption within ±600kms~^ of the center frequency 
(Table 1), to a limit of 30mJy, or 15% of the quasar's miUimeter continuum flux, This 
range includes the z — 0.0523 {cz — 15300 km s~^) redshift seen in absorption against the 
quasar in Na I and Ca II. If the absorbing atomic gas extends 0.3 to 3 kpc on the line of 
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sight, the width of the tidal tail on optical images, then the column densities derived from 
the Na I and Ca II lines by Womble et al. (1990) imply H I densities of 1 - 10cm~^, too 
low for any gas to be in molecular form. 

9.3. 10565+2448 

The main optical and near IR peak of IRAS 10565+2448 has strong Br 7 hues 
(Goldader et al. 1995) and H II region-type hne ratios (Armus, Heckman, & Miley, 1989, 
1990; Veilleux et al. 1995). Murphy et al. (1996) suggest this merger is a triple galaxy 
system. On their r band image there is a secondary source on a tidal tail, 26" (20kpc) 
northeast of the main peak, at the same redshift. A third object with unknown redshift is 
7.9" (6.2 kpc) southeast of the main peak. In our 500 MHz band we detect no CO toward 
the other two r-band sources, to a limit of 20mJy beam~^ in 20kms~^ channels. At 
centimeter wavelengths, there is extended emission and an unresolved, nonthermal core 
(Condon et al. 1991; Crawford et al. 1996). 

Figure 27 shows the maps of CO integrated intensity, isovelocity contours, and CO 
linewidth. The peak of the 2". 3 x 1".7 CO source coincides with the compact cm-radio 
continuum core. Along the line of nodes, the CO line profiles are very symmetric, with 
well-defined blueshifted peaks in the east and redshifted peaks in the west (Fig. 28). 
Together with VII Zw 31 and Arp 193, the galaxy 10565+2448 is one of the best examples 
of a rotating ring in our sample. In the position- velocity diagram (Fig. 29), there are two 
prominent peaks, and a velocity shift of 180kms~^ over the central 1" in R.A. Models that 
best fit the data are those with a minimum ring radius -Rmin in eq-(l) equal to Rq = 230 pc. 
As with VII Zw 31 and Arp 193, filled-disk models, with i?min = 0, give poorer fits, with 
much less contrast. The CO ring must be nearly face-on, because of the small separation 
(80kms~^) of the twin peaks, the narrow CO hnewidth (140kms~^ FWHP), and the small 
extent of the integrated CO intensity on the kinematic major axis. Unlike Arp 220 or 
VII Zw 31, there is not much CO in an outer disk beyond the nuclear ring of 10565+2448. 

9.4. 17208-0014 

The ultraluminous galaxy 17208—0014 has Lir = 2.2 x 10^^ Lq and an optical and 
near-IR line spectrum that indicates H II region-type excitation and high reddening (Martin 
et al. 1989; Kim et al. 1995; Goldader et al. 1995; Veilleux et al. 1995). Images at r-band 
(6550 A) show two tidal tails from a merger (Melnick & Mirabel 1990; Murphy et al. 1996). 
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The innermost parts of the two tails, extending north and east to radii > 20 kpc are nicely 
shown in the r-band image by Sanders & Kim (in Solomon et al. 1997). At 2". 9 (2.2 kpc) 
southeast of the nucleus, the r band and / band (0.82 fim) images show a secondary peak 
on a tidal tail. The i^'-band (2.2 /xm) surface brightness has an r^/^ profile out to a radius 
of 10 kpc, as in elliptical galaxies, probably due to the merger. The nucleus has a size of 
1".8 X 1".4 in K band and appears single, which suggests a completed merger (Zenner & 
Lenzen 1993; Murphy et al. 1996). The nucleus also has one of the strongest known OH 
megamasers, with 10^ in the 1667 MHz OH line, and a radio continuum source of size 
0".32 X 0".26 (220x270pc; Martin et al. 1989). 

Figure 30 shows the maps of CO (1-0) integrated intensity and isovelocity contours. 
The CO source has a size of 1".8 x 1".6, smaller than measured by Planesas, Mirabel, & 
Sanders (1991), and about the same size as the K band source. The centroid of the CO 
source (Table 2) agrees to within 0".2 with the cm-radio continuum source (Martin et al. 
1989). The isovelocity contours and the CO spectra show the molecular gas is blucshiftcd 
in the northwest and redshifted in the southeast. The position-velocity diagram (Fig. 30) 
indicates a strong velocity gradient with a change of 400kms~^ over 1".5 at p. a. 120°. 
We interpret this angle as the line of nodes of the nuclear disk. The CO spectrum at the 
peak of the source (Fig. 30) has a linewidth of 375kms~^ FWHM, and 700kms~^ to zero 
intensity. At the CO peak, there appears to be a weak continuum at the 5 mJy level. 



9.5. 23365+3604 

On i?-band images, the ultraluminous galaxy 23365+3604 has two tidal tails, but only 
a single, blue nucleus (Klaas & Elsasser 1991; Murphy et al. 1996). The optical spectrum 
(e.g., Veilleux et al. 1995) is of the LINER type, with line ratios implying shock velocities 
of 80 - 90kms~^ and pre-shock electron densities of 1 - 10cm~^ in tenuous gas, as in old 
supernova remnants (Klaas & Elsasser 1991). In the nucleus. The [O HI] line has a width 
of 260kms~^ FWHM (Kim et al. 1995). The source has strong Br7 and H2 hues at 2//m 
(Goldader et al. 1995). In the 6 cm radio continuum, there is a secondary source extending 
1" south of the nucleus (Crawford et al. 1996). 

In our 40kms~^ channel maps, the CO source is compact in all channels over a range 
of 320kms~^. The CO source has a size of 1".0 x 0".9, with a peak 2" south of the optical 
position listed by Klaas & Elsasser (1993). The gradient in the isovelocity contours is at p. a. 
135°, and the CO spectra along this line of nodes show asymmetric profiles characteristic 
of a rotating ring or disk. Along this kinematic major axis, the position-velocity diagram 
shows a shift of 300kms~^ over 4" (Fig. 31). 
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lo. DISCUSSION 



10.1. Source sizes, gas masses from CO imaging, and the ratios M^s^/L'qq and 



The main resuhs of this study are the CO source sizes — the halfwidths of the 
integrated CO emission (Table 3) , and the inner and outer disk radii from the model fits to 
the position- velocity diagrams and the channel maps (Table 4). The measured radii of the 
integrated CO correspond to the half-power radii Ri derived from the kinematic data. We 
had previously noted that for ultraluminous galaxies the CO (1-0) flux (Jykms^-*^) was ~ 4 
times the 100 yum flux (Jy), as expected from a blackbody model for the far IR (Downes, 
Solomon, & Radford 1993). We derived the dust temperature from the FIR fluxes, and 
took this to be the intrinsic CO brightness temperature. Prom the CO luminosity, we 
then predicted the mean radius, -Rco, for a spherical CO source (Solomon et al. 1997). 
Interestingly, the sizes measured with the interferometer are within a factor of two of the 
blackbody sizes. The radii Rq and Ri from our fits to the kinematic data (Table 4) bracket 
our single-dish estimate: Rq < Rqo < Ri- The CO radii measured with the interferometer 
now allow us to better estimate the gas and dynamical masses (Table 9). 

In the centers of ultraluminous galaxies, we find M^ss/L'qq 0.8 (Kkms~^ pc^)~^, 
about 5 times lower than the standard value for self-gravitating molecular clouds. We had 
anticipated this result from single-dish data (Downes et al. 1993; Solomon et al. 1997), and 
it is now confirmed by the interferometer maps. We had earlier used the CO linewidths and 
radii from our blackbody model to obtain dynamical mass, which in turn led us to revise 
downward the gas masses from single-dish CO luminosities. We estimated a minimum gas 
mass by assuming optically thin CO. These estimates, which we labeled Mthin (Solomon 
et al. 1997), are within a factor 1.5 of the gas masses derived from our model fits to the 
interferometer data. Table 9 also lists the ratio of gas mass to dynamical mass, which is 
~ 1/6 in the inner, high-density disk, and 1/10 in the outer, low-density disks. 

We showed earlier (Downes et al. 1993) that the Mgas/Mdyn and M^^/L'qq ratios are 
related by 




(10) 



where, letting / = Mgas/Mdyn, we have 



gas 



CO 



(11) 



The factor a = Cn^'^/Ti,, where n is the mean H2 number density over the whole volume 
and Tf) is the CO brightness temperature. In the units used here, the constant C = 2.6 for a 
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sphere and ~ 1.0 for a flared disk with thickness/radius ratio of 0.15, as in some of the disks 
in our sample. These equations depend only on gravity, and are independent of whether the 
gas is self-gravitating or not, or whether the CO lines are optically thick or thin. 

The derived values for Mgas/M^yn and M^ss/L'qq are compatible with starbursts. 
Assume a starburst powers an ultraluminous IR galaxy, with L/M^evr* — 300 LqM©"^ for 
the new stars in the burst, all the power emerging in the IR. Such galaxies typically have 
Lir/L^q = 200Lo(Kkms~^ pc^)~^. Combining the two luminosity ratios gives 

M,ew* = (2/3)L^o > (12) 

or ~ 3 X 10^ Mq as the mass of new stars in the burst. Most galaxies prior to the merger 
would have had a rotation velocity ~250kms^^ at a radius of 500 pc, which implies an old 
stellar bulge about twice as massive as the newly-formed population. Since M^y^ — M^^s + 
-^new* + M3id*; wc havc from the previous equation that 

Mdyn = Mg3. + 2L'co. (13) 
Solving for the molecular gas mass yields 

^ = (l + a2)V2_i . (14) 
-^co 

This equation can be generalized to arbitrary mass ratios of the old and new stellar 
populations. For the Arp 220 disk, with a mean H2 density ~ 500 cm~^ and < T;, > = 20 K, 
these equations predict a ~ 1, M^^/M^y^^ = 0.2, and M^i^/ L'qq = 0.4Mq (Kkms~^ pc^)~^, 
close to the values derived from our model fits to the interferometer data. In other words, 
the starburst scenario is compatible with the derived ratios. 



10.2. Dense gas traced by HCN and extreme steirburst regions 

There is another important gas component that is better traced by emission from high 
dipole moment molecules like HCN, rather than CO. The HCN emission requires densities 
n{B.2)> lO^cm"^. In giant molecular clouds in the disks of the Milky Way and normal 
spirals, the average HCN emission is weak with a typical ratio of Lhcn/-^co= 1/20 to 1/40 
(Gao & Solomon 1998). In the Milky Way, HCN emission is strong only in cloud cores that 
form high-mass stars. We showed previously that ultraluminous galaxies have abnormally 
high HCN luminosities (Solomon, Downes, & Radford 1992), with -Lhcn/-^co= 1/4 to 1/8, 
indicating a higher fraction of the total molecular gas is in dense, > lO^cm"^, star-forming 
cores than in normal galaxies. The densities we derive for the smoothly distributed gas 
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in the rotating disks (see Table 6) are typically 10 times lower than that needed for HCN 

excitation, except for the extreme starburst regions Arp220-west, Arp220-cast and the disk 
of Mrk231. The denser gas responsible for the HCN emission (and CS emission) is thus not 
adequately accounted for in most sources by the "diffuse" CO emitting component that fills 
the disk. 

In Arp 220, the HCN lines we observed have the same velocities as the east and west 
"nuclei", so it seems that the cast and west "nuclei" alone account for most of the HCN 
emission. These dense, compact sources have a hydrogen column density of 0.6 x 10^^ cm~^ 
and mean density of 20,000 cm~^, enough to thcrmalize the lower rotational levels of HCN 
by a combination of collisions and radiative trapping. Within the Arp 220 east and west 
sources, the HCN emission may thus have the same intrinsic brightness temperature as the 
CO (2-1) emission, namely, 50 K (Table 8). Using the sizes and linewidths from Tables 4 
and 5, we estimate the HCN luminosity of these two regions alone to be -Lhcn~ 7 x 10^ 
Kkms"-*^ pc^, which is 3/4 of the observed total (Solomon, Downes, & Radford 1992). These 
two regions thus emit only 1/4 of the CO luminosity but most of the HCN luminosity. 
Their total gas mass of 1.7 xIO^Mq is already accounted for in the mass budget in Table 9. 
Thus for Arp 220, only a small increase in the mass budget would be enough to account for 
the HCN emission. 

While some of the HCN emission undoubtedly comes from dense star forming molecular 
cores embedded in the (relatively) diffuse CO disks of ultraluminous galaxies, it is likely 
that high density, extreme starburst regions similar to Arp 220 east and west exist in most 
of the other ultraluminous galaxies in our sample, and are the real sources of most of the 
HCN emission. Arp 220 is the closest ultraluminous galaxy and has the best resolved 
structure. The second closest galaxy, Arp 193 (see section 6 and Figs. 10 and 11), also 
shows evidence of an extreme starburst in the southeast, that appears similar to the strong 
sources in Arp 220. We suggest that the high-density, extreme starburst regions are the 
source of much of the HCN luminosity but only a fraction of the CO luminosity. Most are 
beyond detection at our current resolution in CO emission since their size of ^100 pc would 
be < 0".14 for all but two of the galaxies (sec Table 3). These objects may then add about 
1 X IO'^Mq, or 25%, to the total gas mass in galaxies with very high HCN emission. This 
is less than we estimated in our 1992 paper. In line with this re-interpretation, we note 
that the galaxy VII Zw 31 has much weaker HCN emission than most of the others in our 
sample. It also has the largest CO disk, a lower ratio of -Lfir/-^co (less star formation per 
solar mass of gas), and may not have any extreme starburst regions like Arp 220 east or 
west. 
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10.3. Stability of the moleculeir disks 

The standard parameter for characterizing the stabihty against local, axisymmetric 
perturbations of a disk that is supported by differential rotation and random motion is 

• (^«) 

where is the one-dimensional random velocity dispersion, k is the local epicyclic 
frequency, and E is the mass surface density. (Safronov 1960; Toomre 1964; Goldreich & 
Lynden-Bell 1965; for a recent review of the criterion applied to gravitationally coupled stars 
and gas in a disk, see Jog 1996). If Q is < 1, the structure is unstable, and large massive 
star clusters may form. We used our interferometer data on the sizes, turbulent velocities, 
rotation velocities, and mass surface densities to estimate the stability of the nuclear disks 
of the ultraluminous galaxies. In our model rotation curves, the epicyclic frequencies are 
K, 2V/R — constant on the rising part of the rotation curve, and k \/2V/R on the flat 
part of the rotation curve. The objects in our sample have high gas mass fractions and 
typically have Qgas < 2, but Qs+g < 1 for stars plus gas (the regime in Fig. le of Jog 1996). 
For Arp 220, we obtain Qgas = 2.2 ai R — 480 pc for the gas alone, but Q ~ 1 for the 
gas plus stars. For VII Zw 31, at a radius of i? = llOOpc, we obtain Q = 1.1 and 0.9 for 
these two masses. These values suggest the molecular disks of Arp 220 and VII Zw 31 are 
globally unstable against axisymmetric perturbations and will form massive star clusters. 

Empirically, this instability in the central disks appears to produce one or more large 
clumps of dense molecular gas — the compact, nearly unresolved peaks on the CO maps. 
Table 12 lists some of these compact regions of dense molecular gas, which we identify with 
extreme starbursts. Their typical radius of 70 to 100 pc may be the scale on which the 
two-component (stars + gas) system becomes unstable. Their mean II2 density can reach 
2 X 10^ cm^^, and their mass appears to be ~ 1 x 10^ Mq of gas initially, which then forms 
~ 10^ OB-type stars over 10^ yr, or ~ 1000 times the number of OB stars in 30 Doradus. 
The total luminosity of one of these extreme starburst regions is 3 x 10^^ to 1 x 10^^ L©. 
The central disk of an ultraluminous IR galaxy typically contains two or three such extreme 
starburst regions at any given time. It is mainly these regions that provide the input 
power to the ultraluminous galaxies. They heat the dust in the central disks to typical 
temperatures of 75 K, (the blackbody fits to the colors measured by IRAS). Since the dust 
is opaque at 100 /im, the disk radiates as a black body. The typical disk radius is 200 to 
300 pc, and the Stefan-Boltzmann formula for the disk as a whole yields a luminosity of the 
order of lO-'^^L©, the typical output power of the ultraluminous galaxies. If these starbursts 
occur in the old, pre-merger nuclei, they change them considerably, creating a new cusp 
with a much greater density of stars. If the extreme starbursts occur slightly outside of the 
old pre-merger nuclei, they will form new nuclei - new cusps of high stellar density. 
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10.4. Why steirs outshine black holes 

The CO disks' radii of ~500pc and rotation speeds of ~300kms^^ yield orbital periods 
of lOMyr. If their FIR luminosity of ~ 10^^ Lq comes from new stars, then our estimates of 
Lpm/Mnew* = 300 LqM©^^ imply star forming rates of SOMoyr^^ The mass of gas plus 
new stars is ~ 6 x 10^ M©, so in 10 rotations of the molecular disks, half the gas turns into 
new stars. 

If part of the ~ 10^^ Lq luminosity came from a black hole accretion disk radiating at 
L ~ 0.1 mc^, the accretion rate would be ~ 1 Mq yr~^, so in 10 molecular disk rotations the 
black hole could only accrete ~ 1% of the molecular gas. The other 99% of the gas would 
continue to form stars over this long period, outshining the black hole. 

Alternatively, if the gas fell in faster and created a 6 x 10^ Mq black hole in one or 
two orbital periods, then the accretion rate would be GOMgyr"^, and a standard Shakura 
& Sunyaev (1973) optically thick accretion disk, radiating at one-tenth the Eddington 
luminosity, would have a luminosity shooting up to 10^^ Lq — 1000 times more luminous 
than the IR ultraluminous galaxies and quasars in the local universe. If the gas fell in 
rapidly at an accretion rate of 60MQyr~^, but the luminosity stayed at 10^^ Lq, then 
the accretion disk's radiative efficiency would be only L = 10~^ mc^ — not any more 
efficient than producing energy from starbursts. Although it is unlikely that such a massive, 
high-density, accretion flow would be advection dominated (e.g., Rees 1982), if the flow 
were in this regime, then almost by definition, the accretion disk would not be a luminous 
object. 

The model of rapid accretion to a black hole thus poses more problems than it solves. 
With a very high accretion rate of GOMQyr"^, why would the IR ultraluminous galaxies 
only have an output of 10^^ Lq if their power source were a black hole? Why aren't they 
as luminous as the powerful quasars at high redshifts, since there is no shortage of fuel? 
Inversely, if the power source were a black hole accreting at a modest rate of IMQyr"^, it 
would last for 6 x 10^ years and there would be many more ultraluminous galaxies than are 
observed. Since the source statistics indicate the ultraluminous phase in mergers lasts for 
~ 10^ years, the more likely answer is that the gas is used up in forming stars - not quasars. 

In summary, even in 10 rotations of the molecular disk, a black hole radiating 10^^ Lq 
at high efficiency could accrete only 1% of the molecular disk's mass. So even though new 
stars convert matter into energy less efficiently than a standard accretion disk, they make 
up for this in ultraluminous IR galaxies by occurring in burst, in a 2000 times larger volume, 
with a 2000 times larger total mass than is available within the Bondi radius GMbh/^^ of 
a supermassive black hole. The molecular disks probably survive for at least 10 rotations. 
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During this time, at most 10" Mq is accreted to a black hole. The rest of the molecular gas 
forms stars. 



11. CONCLUSIONS 

1. ) Rotating disks: At sub-arcsecond resolution, our CO maps of ultraluminous galaxies 
show rotating disks of molecular gas that has been driven into the centers of the mergers. 
The maps of Mrk 273 and Arp 220 also show large-scale streamers or tidal tails roughly 
perpendicular to the nuclear disks. 

2. ) CO only moderately opaque: Model fits show that the observed double-peaked CO 
spectra, and the peak-to-center contrast in the twin-peaked patterns in CO position- velocity 
diagrams, are produced by radiative transfer through the rotating disks, in CO that is only 
moderately opaque. The CO(l-O) opacities are 4 to 10 at the peaks, and lower elsewhere in 
the disks. 

3. ) High turbulence in the molecular disks: In the Arp 220 disk, the 1-D velocity dispersion 
of the molecular gas is a = 100 km s-^ (FWHM 230 km s"^). This high turbulence is one of 
the reasons why the CO is less opaque than in quiescent Milky Way molecular clouds. The 
turbulence also determines the disk thickness, and the heat input to the gas. The turbulent 
heat input, MgasAl^^/i?, is ~ 10^ Lq, the same as the output from the 17;um S(l) line of 
H2, a principal cooling line (Sturm et al. 1996). The line luminosity is that expected from 
IO'^Mq of gas at ~ 100 K. 

4-) Most of the CO luminosity comes from relatively "low-density" gas: Except for 
the compact cores in Mrk 231, Mrk 273, and Arp 220, the molecular disks have true 
CO brightness temperatures ~ 20 K. This is rather lower than the dust and gas kinetic 
temperatures, which are ~65 to 100 K. This implies that the density over most of the 
molecular disk is relatively low (300 to 2000 cm~^), giving subthermal excitation. The 
CO(2-1)/(1-0) ratios are also consistent with the CO being subthermally excited. 

5.) The CO lines in the molecular disks come from a continuous medium, not from 
self- gravitating ( "virialized" ) clouds: At the densities of lO^cm"^ needed to explain the 
observed flux, low-density molecular clouds would be unstable against tidal shear in the 
rotating disks. Because we also detect high density tracers like HCN and CS, we estimate 
~ 10% of the gas is in highly opaque, high-density (10^cm~^), self-gravitating clouds that 
are stable against tidal forces. In the outer disks, however, the CO emission may come from 
"normal" molecular clouds. In this outer disk gas, the volume filling factor would be ~ 0.1, 
the area filling factor 0.3, as in normal galactic disks. 
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6. ) Low dust flux consistent with lower molecular mass: Only in Arp 220 is there significant 
thermal continuum flux from dust at 1.3 mm. The continuum in Mrk 231 at 3 and 1.3 mm 
is nonthermal. At 3 mm, no thermal continuum is detected from dust, to limits of ~2 to 
lOmJy, in accord with the gas masses obtained from CO, corresponding to a CO luminosity 
to gas mass conversion factor about 5 times lower than in self-gravitating clouds in the 
Milky Way 

7. ) Still lots of gas, nevertheless: In spite of the "low" gas densities and CO line opacities, 
the derived gas mass is high; the mass of ~ 5 x 10^ M0 is equal to the mass of molecular 
clouds in a large gas rich spiral galaxy. Within the molecular disks, the ratio of gas 
mass to the enclosed dynamical mass is Mgas/iWdyn = 1/6. The ratio of gas to total 
mass surface density, /i/iitot, reaches a maximum value of 1/3 within the molecular 
disks. The ratio M^ss/L'qq of gas mass to CO luminosity is about one-fifth of its value 
in self-gravitating molecular clouds. For the galaxies in this sample, typical ratios are 
Mg^/L'co = 0.8MQ(Kkms-^pc2)-^ 

8. ) Extreme Starburst Regions: Four extreme starbursts are identified in the 3 closest 
galaxies in the sample including Arp 220, Arp 193 and Mrk 273. They are the most 
prodigious star formation events in the local universe, each representing about 1000 times 
as many OB stars as 30 Doradus. They have a characteristic size of only 100 pc, with about 
10^ Mq of gas and an IR luminosity of ~ 3 x 10^^ Lq from recently formed OB stars. 

Arp 220 has 2 extreme starbursts. The integrated CO and 1.3 mm continuum maps of 
Arp 220 show two compact peaks. The west nucleus, at czisi — 5340 km s~^, is the same 
as the cm-radio and IR west peak, the CO east gas at 5330 km s~^, is associated with the 
cm-radio east peak and its OH megamasers. The CO east gas at 5650 km s~^ is associated 
with the ionized gas in the K-hand cast nucleus. The IR luminosity of the compact peaks 
in Arp 220 can be explained by extreme starbursts in their early phases. The cm-radio 
continuum, the CO intensity, and the 1.3 mm dust flux all suggest the currently most active 
starburst is the west peak. Arp 220 west also shows a complex velocity structure which 
may indicate a bar in formation or a huge molecular outflow. The mass of these compact 
regions is dominated by molecular gas and young stars, not by a bulge population. 

9. ) Regions of dense molecular gas are regions of extreme starbursts: We suggest that 
the HCN emission, which is very strong in most ultraluminous galaxies, originates in the 
high-density, extreme starburst regions similar to Arp 220 east and west. These regions are 
larger than ordinary GMCs, but are filled with molecular gas at a density usually found 
only in small cloud cores. They do not produce most of the CO luminosity, but they do 
emit most of the HCN luminosity. This explains the high HCN luminosity and directly 
relates the HCN emission to star formation, as we suggested previously (Solomon et al. 
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1992). 

10.) Ultraluminous galaxies are powered by starhursts, not AGNs: The CO data show the 
gas in ultraluminous galaxies is in extended disks that cannot intercept all the power of 
central AGNs, if they exist. As a rule of thumb, if you can see the AGN in the UV/visible, 
then it is not heavily absorbed, and cannot be responsible for the far-IR/sub-mm luminosity. 

Furthermore, if the rotating molecular disks were very thin (30 pc) and extended all the way 
in to central AGNs, then their dust would be very hot, emitting most of their power at 10 
to 20 /xm. In fact, most of the dust (and gas) observed in the ultraluminous IR galaxies is 
cool (70 K), emitting at 60 to 100 yUm — the usual temperature of dust in molecular clouds 
heated by starbursts. We conclude that in ultraluminous galaxies — even in Mrk 231 that 
hosts a quasar — the far IR luminosity is powered by extreme starbursts in the molecular 
disks, not by dust-enshrouded quasars. 
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Fig. 1. — VII Zw 31 : left: CO(l-O) in40kms ^ channels. Radial velocities (kms ^, upper 
left of each box) are relative to 109.380GHz {cz\sr = 16147kms~^). Contour unit: llmJy 
bcam-\ with Tb/S = 88K/Jy. Beam = 1".4 x 0".8. The cross at (0,0) is 05'^16°H6.^51, 
79°40'12."8 (J2000). middle: Maps from the model, with same contours and resolution as 
the data, right: Data minus model, with same contour unit. 




Fig. 2.— VII Zw 31 : upper left: CO(1~0) integrated over (-40, +260kms-^). Beam 
= 1"A X 0".8. Contour step 1.2 Jy beam-^ kms'^, with T^/S = 88K/Jy. upper right: 
Velocity contours in steps of lOkms"^ relative to 109.380 GHz {czisr — 16147 km s"^). Labels 
are in kms-^ lower left: C0(2-l) integrated over (-140, +80kms-^). Beam l".l x 0".8. 
Contours: 2 to 16 by 1, in units of 1.6 Jy beam~^ kms~^, with Th/S = 30K/Jy. lower right: 
Velocity contours in steps of lOkms"^ relative to 218.650 GHz {czisr = 16300 km s~^). Labels 
are in kms'^ Center position in aU maps is 05^16™46."51, 79°40'12."5 (J2000). 




Fig. 3. — VII Zw 31 : Position-velocity cuts in declination, with zero position at 
05^16'^46.'^51, 79°40'12."5 (J2000). upper left: CO(l-O) in 1".4 x 0".8 beam. Contour 
step: 6mJy beam~^, with T^/S = 88K/Jy. Velocity is relative to 109.380 GHz {czisr = 
16147kms~^). upper right: Predicted CO(l-O) from the model, with same contours and 
resolution as for the data, lower left: C0(2-l) in l".l x 0".8 beam. Contours: 8 to 90, in 
steps of 4.9 mJy beam~^, with T^/S — 30K/Jy. Velocity is relative to 218.650 GHz (c2;isr = 
16300 km s~^). lower right: Predicted C0(2-l) from the model, with same contours and 
resolution as for the data. 
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Fig. 4. — VII Zw 31 : (a), upper: CO(l-O) spectra vs. radial velocity relative to 
109.380 GHz {cz\sr ~ 16147 km s~^). In the upper left of each box are offsets (arcsec) on 
the kinematic major and minor axes. Beam = 1".4 x 0".8, with T^/S — 88K/Jy. 
(b), lower: C0(2-l) spectra vs. radial velocity relative to 218.650GHz (c^isr = 16300 kms~^). 
In the upper left of each box are offsets (arcsec) on the kinematic major and minor axes. 
Beam I'M x 0".8. The (0,0) position in both diagrams is 05^16'^46.'51, 79°40'12."5 (J2000). 
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Fig. 5— Mrk 231 : (a), left: CO(l-O) in 40kms~^ channels. A 63mJy continuum point 

source has been subtracted at the CO centroid (cross). Radial velocities (kms~^, upper 
left of each box) are relative to 110.602GHz {czisr = 12656 kms~^). Contour unit: 20mJy 
beam-^ Beam = 1".3 x I'M (lower right) with Tb/S = 70K/Jy. 

(b), right: C0(2-l) in 40kms^^ channels. A 36mJy continuum point source has been 
subtracted at the CO centroid (cross). Radial velocities (kms^\ upper left of each box) 
are relative to 221.204GHz {czisr ~ 12650kms~^). Contour unit: 25mJy beam~^. Beam 
= 0".7 X 0".5 (lower right) with T^/S = 75 K/Jy. 
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Fig. 6.— Mrk 231 : 

— upper: CO(l-O) integrated intensity, velocity, linewidth (FWHM), and the 3mm 
continuum. The continuum source has been subtracted from the CO data. Integration 
hmits: (-180, +200kms-^). Beam = 1".3 x l".l. Contours: 

integrated CO: 1 to 18 by 1, in units of 2 Jy beam"-*^ kms"-*^, with Tj^/S — 70K/Jy; 

CO velocity: —60 to +100kms~^, in steps of 10kms~^ relative to 110.602 GHz (c2;isr = 

12656 km s~^). Labels are in kms~^; 

CO linewidth: 50 to 190kms~^ in steps of lOkms""^. Labels are in kms~"^. 
3 mm continuum: 1 to 13 by 1, in units of 4mJy. 

— lower: C0(2-l) integrated intensity, velocity, linewidth (FWHM), and the 1.3mm 
continuum. The continuum source has been subtracted from the CO data. Integration 
limits: (±180 km s"^). Beam = 0".7 x 0".5. Contours: 

integrated CO: 1 to 13 by 1, in units of 5 Jy beam~^ kms~^, with T^/S — 75K/Jy; 

CO velocity: —60 to +70kms~^, in steps of lOkms"^ relative to 221.204GHz (c^isr = 

12650 km s~^). Labels are in kms~^; 

CO linewidth: 50 to 180kms~^ in steps of lOkms"-*^. Labels are in kms"-*^. 
1.3 mm continuum: 1 to 4 by 1, in units of 8mJy. 
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Fig. 7. — Mrk 231 : CO (2-1) position- velocity diagram along the line of nodes (p. a. 
90°). Contour levels are 1 to 19, in units of 20mJy bearn"^, with T^/S — 75K/Jy. Beam 
= 0".7 X 0".5. The (0,0) position is 12'^56"14.^260, 56°52'25."13 (J2000). Velocity is relative 
to 221.204 GHz (c^isr = 12650 km s~^). The continuum source has been subtracted from the 
CO data. 



-53- 



Mrk 331 



oh) 

C0(3-l) 



N 



. 1 I I I u I 1 1 1 I I I , 1 1 1 1 1 I I I I 1 1 1 I J I li I U 1 1 1 I I I ; I L ' I U I I 1 1 1 1 ; I L ' I U 



-io. 



I I I I I I 




I I I I L I I I I I I I I 

^.a|ci' 



I I I I I I I I I I I I I 



J 1 1 1 1 1 1 



■(0. 





ii|iii|iii|ii iili l iliiilii j i L iii L iiilii ii|iiiliii|ii j iliii l iiil 

-(0.g|^),3V --(0,-1)^1 ' ' _ ^-0^3,-0,3) ' _ X--0,'s,-(),3) ' _ ^-0,^,-0,30 ' 



S - 



400 
200 





iili l iliiilii J i L iii L iiilii iiliiiiiiilii J i L iii l iiil 




+S00 -300 km s 



E Major axis (p.a. 90°) W 



Fig. 8.— Mrk 231 : C0(2-l) spectra vs. radial velocity relative to 221.204 GHz {cz]sr = 
12650 km s~^; the velocity scale is reversed to emphasize the symmetry of the profiles). In the 
upper left of each box are offsets (arcsec) on the kinematic major and minor axes. The (0,0) 
position is 12^56°'14.^260, 56°52'25."13 (J2000). Beam = 0".7 x 0".5, with Tb/S = 75K/Jy 
The continuum source has been subtracted from the CO data. 
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Fig. 9. Arp 193 : 

— upper: CO(l-O) maps in 40kms^^ channels. Radial velocities (kms~^, upper left of each 
box) are relative to 112.641GHz {czis,- — 7000 kms~^). Contour unit: —5, 5, 10, then 20 to 
170 by 15mJy beam'^ Beam = 1".3 x 0".9 with n/S = 83K/Jy. 

— lower: C0(2-l) maps in 40kms~^ channels. Radial velocities (kms~^, upper left of 
each box) are relative to 225.282GHz (c^isr = 6994kms~^). Contour interval —10, 10 to 
100 in steps of 10, then 10 to 180 in steps of 20mJy beam~^. Beam = 0".6 x 0".4 with 
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Fig. 10.— Arp 193 : Upper left: CO(l-O) integrated from -245 to +270 km s^^ 
Beam = 1".3 x 0".9. Contours: 0.5, 1, 2, 3, 4, 6, then 10 to 45 by 5Jy beam"^ kms-^ 
with Tfj/S = 83K/Jy. Upper right: Velocity contours in steps of 20kms~^ relative to 
112.641GHz {cz\sj- = 7000 kms~^). Labels are in kms^^. Lower left: C0(2-l) integrated 
from -210 to +270 kms^^ Beam = 0".6 x 0".4. Contour step 2 Jy b eam kms with 
Tb/S = 84K/Jy. Lower right: Velocity contours in steps of 20 km s~^ relative to 225.282 GHz 
(c-2isr = 6994 km s~^). Labels are in kms~^; 
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Fig. 11. — Arp 193 : CO position-velocity diagrams along the line of nodes (p. a. 140°), 
left: CO(l-O), contour unit lOmJy beam-\ with T^/S = 83K/Jy. Beam = 1".3 x 0".9. 

Velocity is relative to relative to 112.641 GHz (c^isr = 7000 km s~^). 

right: C0(2-l), contour unit lOmJy beam^^ with T^/S = 84K/Jy. Beam = 0".6 x 0".4. 
Velocity is relative to relative to 225.282 GHz {cz\sr = 6994 km s~^). The (0,0) position in 
both diagrams is 13^20^35.^315, 34°08'22."20 (J2000). 



58 



NE 



HOD 
100 
mJy 



'(-i.i.'iJs) 
- Arp 163 



- CO(l-O) 




Mm i l l M 1 1 1 L M I L M 1 1 1 1 1 1 1 1 L J 1 1 1 1 1 1 1 1 1 1 1 1 L J, M 1 1 1 L 1 1 1 1 1 1 1 1 II 1 1 i rn 

i.a.-o.ff) ' ..[-0,6,-0,8) ' ..((J.^i.e) ' ' ..Cc,«,-a,9r .^i,3,-o,9r 





-300 □ -l-aOO km 3 



SE Major axis {p. a, 140 ) 



N¥ 



NE 



300 



100 



mJy 




-aOO +aOO km s jjajor axis {p.a. 140°) 



N¥ 



-59- 



Fig. 12. — Arp 193 : upper: CO (1-0) spectra. In each box vertical axis is CO intensity 
and horizontal axis is radial velocity relative to 112.641GHz (czisr = 7000 km s"-*^). In the 
upper left of each box are offsets (arcsec) on the kinematic major and minor axes. Beam 
smoothed to 1".7 x 1".6. 

lower: C0(2-l) spectra. Beam = 0".6 x 0".4. Velocity is relative to relative to 225.282 GHz 
(c^isr = 6994kms-^). The (0,0) position in both diagrams is 13^'20"35.'315, 34°08'22."20 
(J2000). 
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Fig. 13.— Mrk 273 : 

a) CO (1-0) maps in 40kms^^ channels. Radial velocities (kms~^, upper left of each box) 
are relative to 111.076GHz {czisr = 11323 kms~^). Contour unit: —3, +3 to 15 by 5, then 
20 to 80 by lOmJy beam~^, negative contours are dashed, zero contour omitted. Beam 
= 1"A X 1".3 with Tb/S = 56K/Jy. 

b) C0(2-l) maps in 40kms^^ channels. Radial velocities (kms~^, upper left of each box) 
are relative to 222.176GHz (c^isr = 11283 kms~^). Contour interval 10 to 85, in steps of 
5mJy beam~^. Beam — 0".6 with T^/S = 66K/Jy. 
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Fig. 14.— Mrk 273 : CO(l-O), integrated over 250 to 350kms-\ relative to 111.076GHz 
{czisr = 11323kms-^). Beam = 1"A x 1".3 with Tb/S = 56K/Jy. Contour step ImJy 
beani"^, negative contours are dashed, zero contour omitted. 
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Fig. 15.— Mrk 273 : Upper left: CO(l-O) integrated over ±340 km s^^ Beam 
= 1".4 X 1".3. Contours: 0.3 to 0.7 by 0.1, then 1, 1.5, 2, 3, 4, 6, then 8 to 36 by 4Jy 

beam~^ kms~^, with T^/S = 56K/Jy. Upper right: Velocity contours in steps of 20kms~^ 
relative to 111.076GHz (czisr = 11323 kms~^). Labels are in kms~^; Lower left: C0(2-l) 
integrated over ±300 km s-^ Beam = 0".6 x 0".6. Contours: 2, 4, 6, then 8 to 32 in steps 
of 4Jy beam~^ kms~^, with Tb/S = 66K/Jy. Lower right: Velocity contours in steps of 
20kms~^ relative to 222.176 GHz {czigy- = 11283 km s~^). Labels are in kms~^. 
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Fig. 16. — Mrk 273 : CO (1-0) spectra. In each box vertical axis is CO intensity, horizontal 
axis is radial velocity relative to 111.076GHz {cz\sr ~ 11323 kms~^). In the upper left of 
each box are R.A. and Dec. offsets (arcsec); the (0,0) position is 13^44"^42.^01, 55°53'13."0 
(J2000). Beam = 2". 9 x 2".l, with n/S = 16K/Jy. 
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Fig. 17.— Mrk 273 : 

(a) CO (1-0) position- velocity diagram in declination, through the center of Mrk 273. Note at 
+300 kms^^ the second CO complex 4" north of the main CO peak. Contour levels: —2, — 1, 
+1, then 2 to 20 by 2, and 20 to 160 by 10, in units of 1.7 mJy beam-^ with T^/S = 16K/Jy. 
Beam = 2".9 x 2".l. Velocity is relative to 111.076GHz {czi^, = 11323 kms"^). The (0,0) 
position is 13^4^42.^01, 55°5313."0 (J2000). 

(b) CO (2-1) position- velocity diagram along the line of nodes of the nuclear disk (p. a. 90°, 
i.e., the R.A. axis). Contour levels: 3 to 15 by 3, then 20 to 70 by lOmJy beam^-*^, 
with Tb/S = 66K/Jy. Beam = 0".6 x 0".6. Velocity is relative to 222.176 GHz {cz^sr = 
11283 kms-i). The (0,0) position is 13M4°'42.'^01, 55°53'13."5 (J2000). 
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Fig. 18.— Arp 220 : a): CO (2-1) integrated intensity, velocity, linewidth (FWHM), and 
the 1.3mm continuum. CO integration limits: (— 320, +300kms~^). Beam = 0".7 x 0".5 . 
Contours: 

integrated CO: 1 to 10 by 1, in units of 11.4 Jy beam~^ kms~^, with Tb/S = 69K/Jy. 

CO velocity: -225 to +225 kms"^ in steps of 25kms"^ relative to 226.422GHz {czisr = 

5450 km s^^). Labels are in kms~^; 

CO linewidth: 25 to 350kms~^ in steps of 25kms^^. Labels are in kms""*". 
1.3mm continuum: 1 to 12 by 1, in units of 5.7mJy beam~^. 

— h): CO (1-0) integrated intensity, velocity, linewidth (FWHM), and east streamer. For all 
maps, the beam is 1".6 x l".l, with T^/ S — 55 K/Jy, and velocities are relative to 113.228 GHz 
{cz\sr — 5410 km s"-*^). Contours: 

CO integrated from -320 to +300kms-^: 0.1, 0.25, 0.5, then 1 to 12 by 1, in units of 11.4 Jy 

beam~^ kms~^; 

CO velocity: —200 to +260 km s^^, in steps of 20kms^^ Labels arc in kms~^; 

CO linewidth: 25 to 350 km s^^ in steps of 25kms^^. Labels are in kms^^. 

CO east streamer, integrated from +85 to +165kms~^ : 0.025 to 0.2 by 0.025, then 0.25, 

then 0.4 to 2.4 by 0.4, in units of 11.4 Jy beam^^ kms~^. 
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Fig. 19. — Arp 220 : — (^aj CO (1-0) maps in 40 km s ^channels. Radial velocities (km s ^, 
upper left of each box) are relative to 113.228 GHz {cz\sr = 5410kms~^). Contours: 1, 2, 4, 
8, 16, 32, in units of 11.4 mJy beam"!. Beam = 1"6 x I'M with T^/S = 55K/Jy. 
— (b) C0(2-l) maps in 40kms~^ channels. Radial velocities (kms~^, upper left of each 
box) are relative to 226.422 GHz {czisr — 5450 km s~^). The cross marks the position of the 
integrated C0(2-l) west peak (Table 2). Beam = 0"7 x 0".5 with T^/S = 69K/Jy. The 
continuum has been subtracted from the CO maps. Contour unit: 50mJy beam~^. On this 
contour scale, the continuum from the west peak would be two contours. 



-70- 




Fig. 20. — Arp 220 : Radio continuum spectrum. The emission below 30 GHz is nonthermal 
(Norris et al. 1985; Beckhn & Wynn-WiUiams 1987; Baan & Haschick 1987; Baan et al. 

1987; Norris 1988; Condon et al. 1991; Sopp & Alexander 1991 and other references in their 
Table 3). Millimeter and sub- mm data show the increasing flux of dust emission (this paper 
(solid circles); Radford et al. 1991b; ScoviUe et al. 1991, 1997; Woody et al. 1989; Carico 
et al. 1992; Eales et al. 1989; and Rigopoulou et al. 1996. The curves are spectra of the 
form z/^'^+"(l + e^^), normalized to the observed fluxes. The free-free opacity, r, has the 
form (z//z/o)~^'\ where i/q is the turnover frequency. The curve for the dust spectrum is for 
an emissivity index n — 1.5. 
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Fig. 21. — Arp 220 : CO (2-1) spectra. In each box vertical axis is CO intensity, horizontal 
axis is radial velocity relative to 226.422 GHz {cz\sr = 5450 km s""*^). In the upper left of 
each box are offsets (in steps of 0."3) on the kinematic major and minor axes, relative to 
15^^34^57.^24, 23°30'11."2 (J2000). Beam smoothed to = I'M x 0".6 at p.a. 34°, with 
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Fig. 22. — a) Arp 220-disk : CO (2-1) position- velocity cuts along the line of nodes (p. a. 
50°). Labels in the upper left of each box are offsets (arcsec) along the minor axis (northwest 
is positive, southeast is negative). The (0,0) position is 15''34™57.'24, 23°30'11."2 (J2000). 
Contour levels are 1 to 17, in units of SOmJy beam~^, with Tft/S" = 44 K/Jy. Beam smoothed 
to = 1".0 X 0".55 at p.a. 45°. 

h) Arp 220-west: C0(2-l) position-velocity cut through the Arp 220-west, along a possible 
major axis at p.a. 300° (see text). Contour interval 17mJy beam~^, with Ti,/ S = 80K/Jy. 
The source peak is 280 mJy beam~^. 

c) Arp 220-east: C0(2-l) posit ion- velocity cut through Arp 220-east at p.a. 50°. Contour 
interval 17mJy beam"^, with T^/S — 80K/Jy. The source peak is 380mJy beam"^. For 
diagrams b) and c, the restoring beam was 0".6 x 0".5 at p.a. 120°, and the (0,0) position 
is 15''34'"57."225, 23°30'11."5 (J2000) (the western continuum peak). In all these position- 
velocity diagrams, velocity is relative to 226.422 GHz {czisr — 5450 km s"-*^). 
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Fig. 23. — (Plate) Arp 220 : CO contours superposed on a false-color presentation of the 
HST V band image of Arp 220 (Shaya et al. 1994). 

— a) CO(l-O) emission mapped with a beam of 1".6 x I'M, with Ti,/ S — 55 K/Jy, integrated 
in two different velocity ranges, relative to c^isr = 5410 km s~^: — left: redshifted CO(l-O) 
emission integrated from +205 to +325 km s"^; contours 0.75 to 1.75 by 0.25, 2, 2.5, 3, 4.5, 6, 
9, 12, 15, 20, 25, in units of 1.14Jy beam~^ kms~^; — right: blueshiftcd CO(l-O) emission 
integrated from —235 to +5kms~^, with contours 0.75, 2.5, 3.5, 6, 12, 24, 36, 50, 60, 70, in 
units of 1.14 Jy beam~^ kms~^. 

h): CO(l-O) emission in the east streamer and C0(2-l) emission: — left: CO(l-O) emission 
in the east streamer integrated from +85 to +165kms~^; contours and beam as in Fig. 23a^, 
left. — right: CO (2-1) emission integrated from —320 to +300kms~^, relative to c^isr = 
5450kms-^ with contours 20, 30 to 180 by 30, then 240, in units of 0.57 Jy beam"^ kms^^ 
beam = 0".7 x 0".5, with Ti,/S = 69K/Jy. 
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Fig. 24. — Arp 220 : Model of the eastern and western nuclei in the molecular disk. The 
rotation velocity is 330kms~^, and the view is pole-on. The arrows indicate the observed 
radial components along the line of sight, at a disk inclination of 40° from facc-on. Velocities 
are relative to a systemic velocity of cz]sr — 5450 km s~^. Labels are kms~^, values in 
parentheses are czisr- 




Fig. 25.— 00057+4021: upper left: CO(l-O) integrated over (-160, +240 kms"^). Beam 
= 2".2xl".l. Contour step = 0.63 Jybeam-^ km s-\ with = 40 K/Jy. upper right: CO 
velocity contours in steps of lOkms"^. Labels are in kms~^. lower left: CO(l-O) position- 
velocity diagram at p. a. 135°. Contours: 3 to 19 by 1, in units of 3.1 mJy beam"^, with 
Tb/S = 40K/Jy. lower right: CO spectrum in the 2". 2 x I'M beam at the source peak. 
Position offsets are relative to 00^'08"20.'58, +40°37'55."5 (J2000), and velocity scales are 
relative to 110.360GHz {czi,, = 13341 kms'^). 




Fig. 26.— 02483+4302 : upper left: CO(l-O) integrated over (+40, +340 km s"^). Beam 
= 2".0 X 2".0. Contour step = 0.61 Jy beam"^ kms'^, with T^/S = 25 K/Jy. upper right: CO 
isovelocity contours from +140 to +220kms~^, in steps of lOkms"^. Labels are in kms~^. 

lower left: CO(l-O) position-velocity diagram at p. a. 0°. Contours: 3 to 14 by 1, in units of 
3mJy beam~^, with T^/S = 32.5 K/Jy. Beam = 2". 5 x 1".3. lower right: CO spectrum at 
the source peak, in the 2".0 x 2".0 beam, with T^/S = 25 K/Jy. Position offsets are relative to 
02^51"^36."01, +43°15'10."8 (J2000), and velocity scales are relative to 109.700GHz {czi^, = 
15225 kms"^). 
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Fig. 27.— 10565+2448 : CO(l-O) integrated intensity, velocity, and linewidth (FWHM). 
Integration limits: (-100, +220 km s~^). Beam = 2".3 x 1".4 (lower right). Contours: 

integrated CO: 1 to 16 by 1, in units of 1.47 Jy beam~^ kms~^, with T^/S = 32K/Jy; 

CO velocity: +20 to +140 kms"^, in steps of lOkms"^ relative to 110.535GHz (c^^isr = 

12846 km s~^). Labels are in kms~^; 

CO linewidth: 30 to 130kms~^ in steps of lOkms"^. Labels are in kms"-*^. 
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Fig. 28. — 10565+2448 : CO(l-O) spectra. In each box vertical axis is CO intensity and 
horizontal axis is radial velocity relative to 110.535 GHz {czisr — 12846 km s~^). In the upper 
left of each box are offsets (arcsec) on the kinematic major and minor axes, relative to the 
CO centroid listed in Table 2. Beam = 2".3 x 1".4, with T^/S = 32K/Jy. 
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Fig. 29. — 10565+2448 : CO(l-O) position-velocity diagrams along the line of nodes (p. a. 
100°). left: — Observed data: Contour levels are 1 to 10, in units of 15mJy beam~^, with 
Tb/S = 32K/Jy. Beam = 2".3 x 1".4. The (0,0) position is listed in Table 2. Velocity 
is relative to 110.535GHz {czisr = 12846 kms~^). right: — Ring/disk model with the 
parameters listed in the Tables. Contours and resolution are the same as for the observed 
data at the left. 




Fig. 30.— 17208-0014 : upper left: CO(1~0) integrated over (-360, +360 kms"^). Beam 
= 5".l X 1".6. Contour step = 5Jy beam~^ kms~^, with Ti,/S — 12K/Jy. upper right: 
CO velocity contours in steps of 20kms~^. Labels are in kms~^. lower left: CO (1-0) 
position-velocity diagram at p. a. 120°. Contours: 10 to 180, in steps of lOmJy beam"^, 
with Tfo/S" = 12K/Jy. lower right: CO spectrum in the 5".l x 1".6 beam at the source peak 
(Table 2). Position offsets are relative to f 7^23°'2f .^^99, -00°17'00."2 (J2000), and velocity 
scales are relative to 110.535 GHz (c^isr = 12846 kms~^). 




Fig. 31.— 23365+3604 : upper left: CO (1-0) integrated over (-240, +180 km s"^). Beam 
= 3". 9 X 2". 5. Contour step - 2Jy beam"^ kms^"*^, with T^/S = IIK/Jy. upper right: CO 
velocity contours in steps of lOkms"^. Labels are in kms~^. lower left: CO(l-O) position- 
velocity diagram at p. a. 135°. Contour interval: 3mJy beam~^. lower right: CO spectrum 
at the source peak, in the 3". 9 x 2". 5 beam. Position offsets are relative to 23''39™01.'^25, 
36°21'08."4 (J2000), and velocity scales are relative to 108.292GHz {czi^, = 19321 kms"^). 
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Table 1. OBSERVING PARAMETERS 



Son Tpp 


Date 




J. iicxov; 


Flux 


(llppin npj^TTi 


J. i \ Ivv o^^CllC 








(GR?) 


pall nrafor 


(Iv) 


( f^VO^f^O T\ Pi ] 


(K Tv-M 


00057+4021 


06/94- 


-11/94 


110.360 


0003+380 


0.5 


2.2 X 1.1, 47° 


40 


02483+4302 


06/94- 


-11/94 


109.700 


0234+285 


1.4 


2.2 X 1.7, 0° 


27 


VIIZw31 


03/96- 


-05/96 


109.380 


0355+508 


2.6 


1.4 X 0.8, 58° 


88 




03/96- 


-05/96 


218.650 


0355+508 


1.0 


0.7 X 0.4, 46° 


78 


10565+2448 


06/94- 


-11/94 


110.535 


0923+392 


5.5 


2.3 X 1.4, 25° 


32 


Mrk231 


11/95- 


-02/96 


110.602 


0923+392 


4.4 


1.3 X 1.1, 0° 


70 




11/95- 


-02/96 


221.204 


0923+392 


2.7 


0.7 X 0.5, 52° 


75 


Arp 193 


06/96- 


-02/98 


112.641 


1308+326 


1.5 


1.3 X 0.9, 37° 


83 




05/97- 


-02/98 


225.282 


1308+326 


0.8 


0.6 X 0.4, 33° 


84 


Mrk273 


03/96- 


-05/96 


111.076 


1308+326 


1.4 


1.4 X 1.3, 116° 


56 




03/96- 


-05/96 


222.176 


1308+326 


0.7 


0.6 X 0.6, 0° 


66 


Arp 220 


11/95- 


-03/96 


113.228 


3C345 


4.8 


1.6 X 1.1, 25° 


55 




11/95- 


-03/96 


226.422 


3C345 


2.4 


0.7 X 0.5, 28° 


69 


17208-0014 


06/94- 


-10/94 


110.535 


1741-038 


2.6 


5.1 X 1.6, 47° 


12 


23365+3604 


05/96- 


-08/96 


108.292 


3C454.3 


5.0 


3.9 X 2.5, 45° 


11 



Note: The frequency is the zero of our velocity scales, not necessarily the CO line center. 



Table 2. CO POSITIONS AND FLUXES 



Source 


CO 


R.A. 


Dec. 


CO flux 


CO flux 


Continuum 


name 


line 


J2000 


J2000 


interferon!. 


30 m tel. 


freq.: flux 






( h m s) 


( ° ' ") 


(Jy kms"-*^) 


(Jy kms"-*^) 


(GHz): (mJy) 



00057+4021 


(1- 


-0) 


00 


08 


20 


58 


+40 


37 


55 


5 


A '~T 

47 


A 

45 


110: 


< 10 


02483+4302 


(1- 


-0) 


02 


51 


36 


01 


+43 


15 


10 


8 


28 


25 


109: 


< 10 


VIIZw31 


(1- 


-0) 


05 


16 


46 


49 


+79 


40 


12 


6 


87 


92 


109: 


< 2 




(2- 


-1) 


05 


16 


46 


53 


+79 


40 


12 


4 


135 


280 


219: 


< 10 


10565+2448 


(1- 


-0) 


10 


59 


18 


15 


+24 


32 


34 


4 


68 


71 


111: 


< 2 


Mrk231 


(1- 


-0) 


12 


56 


14 


21 


+56 


52 


25 


1 


68 


97 


111: 


63 




(2- 


-1) 


12 


56 


14 


22 


+56 


52 


25 


1 


276 


280 


221: 


36 


Arp 193 


(1- 


-0) 


13 


20 


35 


32 


+34 


08 


22 


2 


161 


162 


110: 


< 5 


disk 


(2- 


-1) 


13 


20 


35 


32 


+34 


08 


22 


2 


360 




228: 


7 


SE core 


(2- 


-1) 


13 


20 


35 


35 


+34 


08 


21 


7 


90 




228: 


3 


sum 


(2- 


-1) 


















450 


415 


228: 


10 


Mrk273 


(1- 


-0) 


13 


44 


42 


12 


+55 


53 


13 


5 


78 


86 


111: 


11 


disk 


(2- 


-1) 


13 


44 


42 


10 


+55 


53 


13 


4 


197 




222: 


<10 


core 


(2- 


-1) 


13 


44 


42 


12 


+55 


53 


13 


4 


34 




222: 


8 


sum 


(2- 


-1) 


















231 




222: 


8 


Arp 220 


(1- 


-0) 


15 


34 


57 


24 


+23 


30 


11 


2 


410 


490 


113: 


41 


west 


(2- 


-1) 


15 


34 


57 


221 


+23 


30 


11 


5 


130 




229: 


90 


east 


(2- 


-1) 


15 


34 


57 


29 


+23 


30 


11 


3 


220 




229: 


30 


disk 


(2- 


-1) 


15 


34 


57 


24 


+23 


30 


11 


2 


750 




229: 


55 


sum 


(2- 


-1) 


















1100 


1100 


229: 


175 


17208-0014 


(1- 


-0) 


17 


23 


21 


95 


-00 


17 


00 


7 


132 


162 


111: 


< 5 


23365+3604 


(1- 


-0 


23 


39 


01 


25 


+36 


21 


08 


4 


30 


47 


108: 


< 5 



Errors: Positions at CO(l-O): ±0".2, at C0(2-l): ±0".l ; Fluxes: ±20%; 
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Table 3. MEASURED DIAMETERS OF THE INTEGRATED CO EMISSION 



Source 


Angular 


P.A. 


Ang.size 


Redshift 


Linear 


Semi- 


Semi- 


name 


diameter 


E of N 


distance 




scale 


maj.axis 


min.axis 




(arcsec) 


(deg.) 


(Mpc) 


(kms~^) 


( Pc/" ) 


(pc) 


(pc) 


00057+4021 


l.lx < 0.6 


-22 


165 


13390 


802 


440 


<240 


02483+4302 


1.8 X 1.7 


-4 


188 


15419 


913 


820 


780 


VIIZw31 


2.4 X 2.1 


-5 


196 


16260 


949 


1140 


1000 


10565+2448 


1.6 X 1.4 


-26 


160 


12923 


776 


620 


540 


Mrk231 


0.9 X 0.8 


77 


157 


12650 


761 


340 


300 


Arp 193-disk 


2.8 X 0.8 


-40 


90 


7000 


436 


610 


170 


Arp 193-SE core 


1.1 X 0.5 


-40 


90 


7000 


436 


240 


100 


Mrk 273-tail 


3.1 X 2.8 


5 


142 


11324 


686 


1060 


960 


Mrk273-disk 


0.9 X 0.6 


18 


142 


11324 


686 


310 


210 


Mrk 273-core 


0.35X <0.2 


76 


142 


11324 


686 


120 


<69 


Arp 220-disk 


2.0 X 1.6 


50 


70 


5450 


341 


340 


270 


Arp 220-wcst 


0.31 X 0.28 


-14 


70 


5340 


341 


53 


48 


Arp 220-east 


0.9 X 0.9 




70 


5650 


341 


150 


150 


17208-0014 


1.8 X 1.6 


7 


159 


12837 


771 


690 


620 


23365+3604 


1.0 X 0.9 


-15 


231 


19330 


1120 


560 


500 



Diameters are FWHM from elliptical Gaussian fits to u, v data. 
Diameters: ±0".2 ; position angles: ±10°. 
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Table 4. DISK SIZES, FROM POSITION-VELOCITY DATA 



Source 


Linear 


Rotation curve 


half-power 


Outer disk 


Disk 


name 


scale 


turnover radius 


radius 


boundary 


thickness 




(pc/") 


Ro (pc) 


Ri (pc) 


-Rmax (pc) 


H{pc) 


00057+4021 


802 


240 


480 


1000 


50 


02483+4302 


913 


270 


730 


1400 


54 


VIIZwSl 


949 


290 


1100 


3300 


58 


10565+2448 


776 


230 


700 


1600 


47 


Mrk231 


761 


75 


460 


1100 


23 


Arp 193-disk 


436 


220 


740 


1300 


65 


Mrk273-disk 


686 


70 


400 


1900 


42 


Arp220-disk 


341 


200 


480 


1400 


80 


Arp 220-west 


341 


35 


68 




68 


Arp220-east 


341 


50 


110 




80 


17208-0014 


771 


310 


540 


1500 


62 


23365+3604 


1120 


340 


670 


2200 


67 



Disk thickness H is FWHM perpendicular to the equator. 
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Table 5. ROTATION CURVES, FROM POSITION-VELOCITY DATA 





Data: 


Data: 


Model: 


Model: 


Model: 




line of nodes 


apparent 


inclination 


rotation 


turbulent 


Source 


E of N 


rotation 


O°=face-on 


velocity 


velocity 




IJ.CXi. 


''rot S'iJ^-'- ' 


V 


V ^ 

I' rot 


AV 




(deg) 


(kms~^) 


(deg) 


(kms~^) 


(kms~^) 


00057+4021 


135 


75 


20 


250 r'^ 


100r-°-5 


02483+4302 





70 


15 


270 r'^ 


40^0.0 


VIIZw31 





100 


20 


290 r'^ 


30r°° 


10565+2448 


100 


75 


20 


220 


40^0.0 


Mrk231 


90 


60 


10 


345 r'^ 


go ^-0.3 


Arp 193-disk 


140 


175 


50 


230 


40^0.0 


Mrk273-disk 


90 


200 


45 


280 


140^0.0 


Arp220-disk 


40 


212 


40 


330 


140r-o-2 


Arp 220-west 


300 


103 


20 


300 


100r-°-3 


Arp 220-east 


50 


300 


60 


350 


100r-°-3 


17208-0014 


120 


130 


30 


260 


150r-0-2 


23365+3604 


135 


130 


30 


260 


100r-°-2 



Model rotation curve: = 1 for r < 1 and = for r > 1, and r = R/Rq, 
where Rq is the rotation curve turnover radius (Table 4). 

Turbulent velocity = halfwidth to (1/e) level. AV = 0.6 FWHM = 1.4 a, 
where a = 1-D velocity dispersion. 

Turbulent velocities are valid for Rq < R < i?max (Table 4). 



Table 6. MODEL TEMPERATURES AND DENSITIES 



Source 


Gas kinetic 


< n(H2) > 


< n(H2) > 


name 


temperature 


at i?o 


outer disk 




(K) 


(cm"^) 


(cm"^) 


00057+4021 


72 r-°-5 


400 


40 


02483+4302 


50 r-°-^ 


250 


20 


VIIZw31 


50 r-°-5 


450 


10 


10565+2448 


60 


1200 


30 


Mrk231 


100 r-o-^ 


3600 


60 


Arp 193-disk 


40 


550 


10 


Mrk273-disk 


70 r-°-5 


1800 


60 


Arp220-disk 


64 r-°-5 


900 


30 


Arp 220-west 


150 r-°-5 


22000 


900 


Arp220-east 


100 r-°-^ 


20000 


900 


17208-0014 


70 r-°-5 


600 


60 


23365+3604 


67 r-°-5 


200 


20 



In the scaling laws, r = R/Rq where Rq is the turnover 
radius of the rotation curve (Table 4). At large radii, we 
adopted minimum gas kinetic temperatures of 10 K. 

Gas densities < n{B.2) > are averaged over the inner and 
outer disk volumes. 



-91 - 



Table 7. CO (1-0) LINE PARAMETERS 





Model: 


Model: 




Model: 


Data: 


Source 


CO(l-O) 


True Tb 


Beam 


Predicted 


Observed 


name 


opacity 


at Rq 


FWHM 


Tf, in beam 


Tb in beam 




at Rq 


(K) 


(arcsec) 


(K) 


(K) 


00057+4021 


4.8 


16 


2.2 X 1.1 


2.7 


2.4 


02483+4302 


9.7 


11 


2.2 X 1.7 


1.9 


1.4 


VIIZw31 


7.8 


18 


1.4 X 0.8 


7.2 


7.4 


10565+2448 


9.8 


36 


2.3 X 1.4 


4.6 


4.5 


Mrk231 


3.4 


56 


1.3 X 1.1 


13.5 


14.0 


Arp 193-disk 


9.9 


25 


1.7 X 1.6 


5.0 


6.9 


Mrk273-disk 


3.7 


36 


1.4 X 1.3 


4.8 


4.5 


Arp220-disk 


5.9 


29 


1.6 X 1.1 


10.2 


10.2 


Arp 220-west 




50 


1.6 X 1.1 


10 




Arp 220-east 




50 


1.6 X 1.1 


10 




17208-0014 


3.9 


21 


5.1 X 1.6 


2.5 


2.2 


23365+3604 


6.0 


10 


3.9 X 2.5 


1.0 


0.8 



Model CO abundance: Xco/ (dv/dr) = 8 x 10"^ 
All brightness temperatures are Ray leigh- Jeans. 
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Table 8. CO (2-1) LINE PARAMETERS 





Model: 


Model: 




Model: 


Data: 


Data: 


Source 


C0(2-l) 


True Tfe 


Beam 


Predicted 


Observed 


C0(2-l) 


name 


opacity 


at Rq 


FWHM 


Tf, in beam 


Tft in beam 


to (1-0) 




at Rq 


(K) 


(arcsec) 


(K) 


(K) 


ratio in 3 mm beam 


VIIZw31 


25 


12 


1.1 X 0.8 


5.4 


5.4 


0.7 ±0.2 


Mrk231 


12 


51 


0.7 X 0.5 


30.3 


28.5 


1.0 ±0.2 


Arp 193-disk 


54 


13 


0.7 X 0.6 


8.2 


8.1 


0.5 ± 0.2 


Mrk273-disk 


11 


34 


0.6 X 0.6 


14 


6.6 




Arp 220-disk 


15 


23 


0.7 X 0.5 


16.6 


10.0 


0.6 ±0.2 


Arp 220-west 


17 


53 


0.7 X 0.5 


18.8 


17.0 




Arp 220-east 


9 


70 


0.7 X 0.5 


9.5 


9.1 





Model CO abundance, Xco/ (dv/dr) = 8 x 10~^. 
All brightness temperatures are Rayleigh- Jeans. 
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Table 9. CO (1-0) LUMINOSITY, GAS MASS, AND DYNAMICAL MASS 



Source 


radius 


f 

CO 


-''-'gas 


gas / -^CU 




Mass ratio 


U 1 Lifr^f 


name 


pc 




lO'^ M0 


Mq/L^ 


lO'^ M0 




max. 


00057+40 to Ri 


480 


2.3 


0.6 


0.3 


7.0 


0.09 


0.31 


to Rmax 


1000 


4.0 


1.4 


0.4 


14.5 


0.10 


— 


Uz48c)+4o to Ki 


Ton 


n 
z.U 


i.z 


U.D 


iz.4 


U.iU 




to -Rmax 


1400 


3.2 


1.9 


0.6 


23.6 


0.08 





VIIZw31 to Ri 


1100 


3.2 


3.0 


0.9 


21.5 


0.14 


0.22 


to K2 


1 nnn 

lyuu 


0.0 


o.z 


i.U 


O/.O 


U.i4 




to -Rmax 


3670 


8.5 


8.3 


1.0 


71.5 


0.12 


— 


10565+24 to Ri 


700 


2.3 


1.8 


0.8 


7.8 


0.23 


0.28 


to K2 


1 1 ttn 
iiDU 


A 
0.4 


o.U 


u.y 


io.i 


U.zo 




to -Rmax 


1600 


5.5 


4.0 


0.7 


17.4 


0.23 





Mrk231 to i?i 


460 


2.8 


1.8 


0.7 


12.7 


0.14 


0.20 


to -fl2 


ooU 


4.Z 


Q 1 
O.i 


U. / 


ZO.i 


U.io 




to -Rmax 


1700 


5.1 


4.0 


0.8 


31.5 


0.13 


— 


Arp 193 to i?i 


740 


1.4 


1.7 


0.8 


9.1 


0.19 


— 


to R2 


1100 


2.1 


2.6 


0.8 


15.0 


0.16 


— 


to -Rmax 


1300 


2.6 


3.4 


0.9 


16.1 


0.21 


— 


Mrk273 to i?i 


400 


4.1 


2.1 


0.5 


8.7 


0.24 


0.28 


to i?2 


900 


4.6 


2.7 


0.6 


16.2 


0.17 




to -Rmax 


1900 


5.6 


5.6 


1.0 


34.6 


0.16 




Arp 220 west 


68 


0.6 


0.6 


1.0 


1.4 


0.43 




Arp 220 east 


150 


0.8 


1.1 


1.3 


3.2 


0.35 




Arp 220 disk to i?i 


480 


2.7 


2.0 


0.8 


12.1 


0.17 


0.30 


Arp 220 disk to R2 


760 


3.7 


3.0 


0.8 


19.2 


0.16 




Arp 220 total 


1360 


5.9 


5.2 


0.9 


34.5 


0.15 
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Table 9 — Continued 



Source 


radius 


T' 

^CO 






Mdyn 


Mass ratio 


/^//^tot 


name 


pc 




10^ Mq 




10^ Mq 




max. 


17208-00 to Ri 


540 


3.3 


2.4 


0.7 


8.5 


0.28 




to Rraax 


1500 


5.8 


6.1 


1.0 


23.5 


0.26 




23365+36 to Ri 


670 


3.0 


1.5 


0.5 


10.5 


0.14 




to Rmax 


2200 


5.4 


3.8 


0.7 


34.5 


0.11 





L; = K kms ^ pc^; [i —gas (H2+ He) surface density in MqPC ^; //tot = gas + stars; /////tot 
= peak ratio in the disk. 

Ri = -Rmin + AJ? is the inner disk's half-power radius; R2 = Rmin + 2Ai? is its zero-power 
radius; i?max is the boundary of the outer disk (eq. 1; Table 4). 
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Table 10. CONTINUUM EMISSION 



2.6 mm 2.6 mm 1.3 mm 1.3 mm Interpretation 

Source Predicted Observed Predicted Observed of observed 

name dust flux total flux dust flux total flux flux at 

(mJy) (mJy) (™Jy) (i^Jy) 1.3 mm 



00057+4021 


0.9 


< 10 


02483+4302 


0.6 


< 10 


VIIZw31 


0.9 


< 2 


10565+2448 


2.2 


< 2 


Mrk231 


2.4 


63 


Arp 193 


3.8 


< 5 


Mrk 273-core 


1.5 


11 


Arp 220-west 


6.3 




Arp 220-east 


2.3 




Arp220^isk 


4.8 




Arp220-total 


13.4 


41 


17208-0014 


3.5 


5 


23365+3604 


1.0 


< 5 



10.1 — — 

6.1 — — 

17.0 < 20 — 

24.9 — — 

7.0 36 80% nt, 20% d 

39 10 >90% d 

16 8 50% nt, 50% d 

71 90 10% nt, 10% ff, 80% d 

25 30 10% nt, 10% ff, 80% d 

54 55 10% nt, 90% d 

150 175 15%nt, 5%ff, 80% d 



40 
11 



Predicted fluxes are for the temperatures in Table 6 and dust emissivity index 1.5 
Errors in observed fluxes: ±20%; nt = nonthermal, ff = free-free, d = dust. 
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Table 11. STARBURST CONSISTENCY CHECK: MASS IN GAS, NEW STARS, AND 
OLD BULGE STARS IN THE CENTRAL MOLECULAR DISKS 





Mrk 231 
disk 


Mrk 273 
disk 


Arp 220 
disk 


Reference radius: 
Ripe) 


460 


400 


480 


Luminosity: 
Lfir(< R) (10^2 


1.9 

500 


1.0 

500 


0.7 
300 


Deduced mass in new stars: 

Mnew*(< R) (10^ Mq) 


3.8 


2.0 


2.3 


Mass from CO model: 
Mdyn(< R) (10^ Mo) 
Mga.(< R) (109 Mo) 


12.7 
1.8 


8.7 
2.1 


12.1 
1.9 


Deduced mass in old stars: 
Moid.(< i?) (10^ Me) 


7.1 


4.6 


7.9 



Lia/Mnew* = assumed luminosity-to-mass ratio in a rapid 
starburst. 

= Lfir(< R) / (Lin/Mnewk) = mass of new stars formed 
in the current starburst. 

-^oid* — -^dyn " -^gas ~ -^new* — Hiass of bulge stars withiu 
R{pc) 
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Table 12. PROPERTIES OF EXTREME STARBURST REGIONS 



Arpl93 Mrk273 Arp220 Arp220 
SE core core west east 



1) Reference radius: 
i?(pc) 

2) Gas Mass: 
Mga.(< R) (109 Mo) 

3) Mean gas density: 
< N{R2) > (cm-3) 



150 



0.6 



4) Total Mass (Gas plus stars): 
Mtot(<i?) (10^ Mo) 1.4 



120 



1.0 



68 



0.6 



2.6 



1.4 



110 



1.1 



2 X 10^ 5 X 10^ 2 X 10^ 8 x 10^ 



3.2 



5) Estimated mass in new stars: 

Mnew*(< R) (10^ Mo) 

6) Luminosity: 
LFm{<R) (10^2 Lo) 

7) Luminosity to mass ratio: 
Lpin/M^,^^ (Lo/ Mo) 



0.8 



0.2 



300 



1.6 



0.6 



360 



0.8 



0.3 



380 



2.1 



0.2 



100 



1 ) Radius from measurements (Table 3 or 4) . 

2) Mgas a; I.OL'qq (Table 9). 

3) < N{R2) >{cm-' = 0.1(Mga./Mo) r^, (cylinder). 



4) Mtot = RAvyc. 



Mtot - M,. 



6) Lir/Lo = 0.13r;^r^ with T from Table 6. 

7) (ivprn/Mnew*) = luminosity-to-mass ratio of new stars formed in the 
current star burst. 



This figure "fl8a.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/9806377vl 



This figure "fl8b.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/9806377vl 



This figure "f23a.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/9806377vl 



This figure "f23b.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/9806377vl 



